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METHOD FOR DETECTING HEAT EXCHANGER 
TUBE FAILURES AND THEIR LOCATION WHEN USING INPUT/LOSS 
PERFORMANCE MONITORING OF A RECOVERY BOILER 

[0001] This application is a Continuation-In-Part of pending U.S. Patent Application Serial 
Number 10/268,466 filed October 9, 2002, which issued on November 18, 2003 as U.S. Patent 
6,651,035 for which priority is claimed and is incorporated herein by reference in its entirety. 
Application Serial Number 10/268,466 which, in turn, is a Continuation-In-Part of pending U.S. 
Patent Application Serial Number 10/131,932 filed April 24, 2002, for which priority is claimed 
and is incorporated herein by reference in its entirety. Application Serial Number 10/131,932 
which, in turn, is a Continuation-In-Part of U.S. Patent Application Serial Number 09/273,71 1 filed 
March 22, 1999, which issued on February 18, 2003 as U.S. Patent 6,522,994 and is incorporated 
herein by reference in its entirety; Application Serial Number 09/273,711 which, in turn, is a 
Continuation-In-Part of U.S. Patent Application Serial Number 09/047,198 filed March 24, 1998 
now abandoned. 

[0002] Application Serial Number 10/13 1,932 is also a Continuation-In-Part of U.S. Patent 

Application Serial Number 09/630,853 filed August 2, 2000, which issued on June 24, 2003 as 
U.S. Patent 6,584,429 and is incorporated herein by reference in its entirety; Application Serial 
Number 09/630,853 claims the benefit of U.S. Provisional Patent Application Serial Number 
60/147,717 filed August 6, 1999. 

[0003] Application Serial Number 10/131,932 is also a Continuation-In-Part of pending 

U.S. Patent Application Serial Number 10/087,879 filed March 1, 2002 and is incorporated herein 
by reference in its entirety; Application Serial Number 10/087,879 which, in turn, is a 
Continuation-In-Part of U.S. Patent Application Serial Number 09/273,711 now U.S. Patent 
6,522,994; Application Serial Number 09/273,71 1 which, in turn, is a Continuation-In-Part of U.S. 
Patent Application Serial Number 09/047,198 now abandoned. Application 10/087,879 is also a 
Continuation-In-Part of U.S. Patent Application Serial Number 09/630,853 now U.S. Patent 
6,584,429; Application Serial Number 09/630,853 claims the benefit of U.S. Provisional Patent 
Application Serial Number 60/147,717. Application 10/087,879 is also a Continuation-In-Part of 
U.S. Patent Application Serial Number 09/827,956 filed April 4, 2001, which issued on May 6, 
2003 as U.S. Patent 6,560,563 and is incorporated herein by reference in its entirety; Application 
Serial Number 09/827,956 which, in turn, is a Continuation-In-Part of U.S. Patent Application 
Serial Number 09/759,061 filed January 11, 2001 now abandoned; Application Serial Number 
09/759,061 which, in turn, is a Continuation-In-Part of US Patent Application Number 09/273,71 1 
now U.S. Patent 6,522,994; Application Serial Number 09/273,711 which, in turn, is a 



Continuation-In-Part of U.S. Patent Application Serial Number 09/047,198 now abandoned. 
Application 10/087,879 is also a Continuation-In-Part of pending U.S. Patent Application Serial 
Number 09/971,527 filed October 5, 2001; Application Serial Number 09/971,527 which, in turn, 
is a Continuation-In-Part of U.S. Patent Application Serial Number 09/273,711 now U.S. Patent 
6,522,994; Application Serial Number 09/273,71 1 which, in turn, is a Contmuation-In-Part of U.S. 
Patent Application Serial Number 09/047,198 now abandoned; Application Serial Number 
09/971,527 is also a Continuation-In-Part of U.S. Patent Application Serial Number 09/630,853 
now U.S. Patent 6,584,429; Application Serial Number 09/971,527 is also a Continuation-In-Part 
of U.S. Patent Application Serial Number 09/827,956 now U.S. Patent 6,560,563; Application 
Serial Number 09/827,956 which, in turn, is a Continuation-In-Part of U.S. Patent AppUcation 
Serial Number 09/759,061 now abandoned; Application Serial Number 09/759,061 which, in turn, 
is a Continuation-In-Part of US Patent Application Number 09/273,71 1 now U.S. Patent 6,522,994; 
Application Serial Number 09/273,711 which, in turn, is a Continuation-In-Part of U.S. Patent 
Application Serial Number 09/047,198 now abandoned. 

[0004] This invention relates to a recovery boiler as used by the pulp and paper industry 

burning black liquor, and, more particularly, to a method for rapid detection of tube failures and the 
location of the affect heat exchanger within the recovery boiler, without need for direct 
instrumentation, thereby preventing serious equipment damage, preventing boiler explosion, 
preventing injury to operators and minimizing repair time on the effected heat exchanger. 

BACKGROUND OF THE INVENTION 

[0005] Although especially applicable to "The Input/Loss Method" as installed at recovery 

boilers burning black liquor, this invention may also be applied to any other of the "Input/Loss 
methods" installed at any thermal system burning a fossil fuel. The teachings of this invention may 
be implemented for monitoring of any thermal system burning a fossil-fiiel, or a thermal system 
burning a mix of fossil fuels and morganic fuels. Such monitoring is assumed to be conducted in a 
continuous manner (i.e., on-line), processing one monitoring cycle after another, each cycle 
includes determining stoichiometric balances of the combustion process and, specifically, the fuel's 
chemistry, heating value, boiler efficiency, system efficiency and evaluation for possible tube 
failures. Specifically, The Input/Loss Method and its associated technologies are described in the 
following U.S. Patents and patent applications (cited above): U.S. Patent 6,522,994 (hereinafter 
termed '994), U.S. Patent 6,584,429 (hereinafter termed '429), U.S. Patent 6,560,563 (hereinafter 
termed '563), U.S. Application 10/087,879 (hereinafter termed '879), U.S. Patent 6,651,035 
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(hereinafter termed '035) and U.S. Application 10/131,932 (hereinafter termed '932); and in any of 
their related provisional patent applications and Continuation-In-Parts. One of the Input/Loss 
methods, a rudimentary method, is described in U.S. Patent 5,367,470 issued November 22, 1994 
(hereinafter termed '470), and in Patent 5,790,420 issued August 4, 1998 (hereinafter termed '420). 
[0006] Conventional steam generators and recovery boilers having gas-to-working fluid 

heat exchangers, may be prone to tube leaks of their working fluid (typically water as liquid or 
steam). These tube leaks represent a potential for serious physical damage to heat exchangers due 
to pipe whip (i.e., mechanical movement) and/or steam cutting of metal given high leakages 
flowing at critical velocities. In some recovery boilers, pressures of the working fluid may exceed 
2300 psia. Given failure of a heat exchanger tube, such fluid will experience many times critical 
pressure ratio as it expands into the combustion gases; that is, mixing with the products of 
combustion at essentially atmospheric pressure. When undetected, the damage from such tube 
failures may range from $2 to $10 million/leak forcing the system down for major repairs. If 
detected early, tube failures may be repaired before catastrophic damage, such repairs lasting only 
several days and costing a fraction of the cost associated with late detection and catastrophic 
damage. Repair times may be fiuther reduced if the location of the heat exchanger which has the 
leak is identified before repairs are initiated. 

[0007] However, an unique situation found with recovery boilers is associated with the pulp 

producing process to which they are integrated. This process involves first de-barking and 
chipping wood; then digesting the wood in an aqueous solution of NaOH and Na2S04 (or other 
sodium-based compounds), forming a "white liquor^'; then heating the brew; then separating the 
pulp from the spent liquor, the spent liquor is termed "black liquor" which consists of organics, 
water and inorganics (mostly sodium); and then the black liquor is burned in a recovery boiler. 
The essential fimction of a recovery boiler is the reduction in the fimiace of sodium sulfate 
(Na2S04, present in the black liquor) to sodium sulfide (Na2S). The eflSciency of this sulfiir 
reduction process is gauged by a "Reduction Efficiency" parameter. Heat from combustion of the 
fired organics, originating from the wood digesting process, generates steam. Black liquor 
inorganics, after reduction, are collected at the bottom of the fiimace as a molten smelt, removed 
and recycled to recover sodium. Given a high Reduction EfiBciency, smelt principally consists of 
Na2C03, Na2S, inerts and free carbon. 

[0008] The problem of tube failures in recovery boilers, in addition conventional problems 

cited above, is when water comes in contact with the molten smelt (typically at over 1400 F, with 
heavy concentrations of sodium); explosion is likely and may occur within minutes after tube 
failure. Recovery boiler explosions have dogged the pulp and paper industry since inception of the 
pulp producing process (i.e., called the Kraft process). Recovery boiler explosions injure and kill 
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people every year. From 1948 through 1990 the industry recorded 140 recovery boiler explosions, 
three-quarters of which were smelt-water explosions. To place emphasis on the problem, the 
industry ranks explosions by severity: by definition just a "moderate explosion" keeps the plant 
off-line from 10 to 50 days; whereas a severe explosion keeps the plant off-line more than 50 days 
(typically lasting more than 120 days). 

[0009] As common with conventional steam generators, tube failures in recovery boilers 

are typically caused by one the following general categories: 
Weld failure of heat exchanger tubes; 

Metallurgical damage caused by hydrogen absorption in the metal resulting 
in either embrittlement or the formation of non-protective magnetite; 

Caustic gouging caused by the presence of free hydroxide in the water; 

Corrosion-fatigue damage from the water-side of the tube, compounded by stress; 

Corrosion damage caused by impacts from solid ash particles; 

Fatigue failure caused by oxidation and/or mechanical movement, 
compounded by stress; 

Overheating (e.g., from tube blockage) causing local creep; and 

Physical damage from steam cutting and/or mechanical movement 
associated with another failed tube in the same locale. 
Commonly, the physical leak initiates as a relatively small penetration, although initial breaks may 
also occur. For reference and further discussion see: Chapter 18, "Failure Analysis and In-Service 
Experience - Fossil Boilers and Other Heat Transfer Surfaces" of The ASME Handbook on Water 
Technology for Thermal Power Systems, P. Cohen, Editor, The American Society of Mechanical 
Engineers, New York, NY, 1989; and J. Gommi, "Root Causes of Recovery Boiler Leaks", 1997 
Engineering and Papermakers Conference, TAPPI Proceedings, available from TAPPI Press as 
product code ENG97509, Atlanta, Georgia. 

[0010] Present industrial art associated with conventional steam generators and recovery 

boilers have practiced the detection of tube failures using one or more of six general methods: 1) 
operator interface; 2) acoustic monitoring; 3) water balance testing; 4) monitoring of effluent 
moisture using instrumentation located at the system's effluent boundary (i.e.. Stack); 5) 
monitoring the concentration of chemicals added to the working fluid whose change is sensitive to 
leakage; and 6) through use of artificial neural network technologies. Operator interface involves 
the use of his/her knowledge, experience, listening skills, and visual skills using remote cameras. 
However, all operators do not have the same high skill-set required. Acoustic devices detect the 
unique noise created by fluids at high velocities. However, acoustic devices rarely work in large 
steam generators, are expensive and require benchmarking with known acoustical signatures. 
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Water balance testing may be conducted periodically on the entire system through which large 
water losses due to tube failures might be discovered. However, water balance testing is expensive, 
insensitive to small leaks, and typically may not be conducted at suflBcient frequency to prevent 
serious damage. The use of an effluent moisture instrument has been shown to be sensitive to tube 
failures. Effluent moisture instrumentation may not differentiate between originating sources of 
water (e.g., between high humidity in the combustion air, or high fiiel water, or changing fiiel 
water, or a tube leakage). However, it might be practical to detect tube leakage by monitoring the 
difference in signals, or the rate of change of the difference in signals, between an effluent moisture 
instrument and one monitoring ambient air. Note that typical black liquor fuels contain up to 35% 
fuel water (approximately the same amount of water as found in some Powder River Basin (PRB) 
coals); thus changes in effluent moisture, even referenced to an ambient measurement, may be 
insensitive to small tube failures. 

[0011] Monitoring the concentration of chemicals added to the working fluid operates by 

making a chemical mass balance on the working fluid based on a combination of flow 
measurements and chemical concentration measurements. Computed is a mass balance of a 
specific stable and non- volatile species (such as phosphate or molybdate) which has been uniquely 
added to the working fluid of the boiler. Basically a foreign chemical is injected into the working 
fluid; when a tube leak occurs the concentration of the chemical will change, thus detection. This 
method, developed by Burgmayer, Hong and Gunther, is described in U.S. Patent 6,484,108 issued 
November 19, 2002. A similar method is described in U.S. Patent 3,522,008 issued July 28, 1970. 
A similar method is also described in U.S. Patent 5,320,967 issued June 14, 1994. None of these 
methods involve combustion gases nor any stoichiometric balance involving the combustion 
process. The most serious limitation of these methods is their lack of sensitivity. Burgmayer' s. 
Table 4 presents results of actual tests indicating that for a boiler producing 500,000 Ib/hr of steam, 
detection of a 0.76% leakage (3,800 Ib/hr) took 45 minutes, while a 2.8% leakage (14,000 Ib/hr) 
took 15 min. to detect. Such sensitivities are not adequate to safeguard operators and equipment. 
[00121 U.S. Patent 6,192,352 by Alouani issued February 20, 2001 (and his US Patent 

Application Publication 2001/0001 149 of May 10, 2001) discloses a method to detect tube failures 
using artificial neural network and fiizzy logic technology (ANN). No where in Alouani's patent is 
explicit thermodynamic modeling taught. Alouani's patent teaches that ANN technology may 
learn to predict tube failures through recognition of patterns in raw data. Such raw data may 
include coal pulverizer flow (fiiel flow), boiler drum pressure, reheat temperature, burner tilt 
positions, etc. The disadvantage to this method is that it requires a database fi'om which it may 
leam. Such a database, as must be associated with an actual tube leak, is not fi^uent, 
inconvenient and could not contain a defined tube leakage flow rate (explicit system mass balances 



are not mentioned nor inferred by Alouani). In addition, Alouani's patent Figure 5 indicates that a 
number of days is required for his system to detect a tube leak. In a survey of critical tube leaks in 
recovery boilers, it was found that approximately half of those leaks for which the time between 
leak initiation and explosion was known, the explosion occurred within 15 minutes after leak 
initiation; approximately 75% of the recorded explosions occurred within 30 minutes. This 
survey's reference is: D.G. Bauer and W.B.A. Sharp, "The Inspection of Recovery Boilers to 
Detect Factors That Cause Critical Leaks", TAPPI Joumal, September 1991, TAPPI Press, Atlanta, 
Georgia. Although there are more than a half-dozen vendors offering one or more of the six tube 
leak detection methods, in practice all known methods suffer serious short-comings and are not 
reliable in detecting early tube failures. 

[0013] The patents '470 and '420 make no mention of heat exchanger tube failures nor their 

detection. Although the technologies of Patents '994, '429 and '563 generally support this 
invention, they make no mention of tube failure detection nor their location. Applications '879, 
'035 and '932 support this invention directly. Although the methods of '879, '035 and '932 are 
usefiil, the present invention further improves these methods and applies them to recovery boilers. 
There is no established art directly related to this invention; there is clear need for early detection of 
tube failures and to determine their location within the recovery boiler system. 

SUMMARY OF THE INVENTION 

[0014] This invention relates to a recovery boiler as used by the pulp and paper industry 

burning black liquor, and, more particularly, to a method for rapid detection of tube failures and the 
location of the affect heat exchanger within the recovery boiler, without need for direct 
instrumentation. Tube failures are detected through use of combustion stoichiometrics, in 
combination with an ability to correct effluent data through use of optimization procedures. The 
location of the failure within a recovery boiler is determined through energy balances, high 
accuracy boiler efficiency and iterative techniques. Further, this invention teaches how the 
stoichiometric mechanism of tube failure may be identified and reported to the system operator. 
This invention addresses the deficiencies found in all present detection methods. 
[0015] Effluent water concentration (at the Stack) may consist of any one or all of the 

following sources of working fluid (assuming the working fluid is water): heat exchanger tube 
leaks; water added at the point of combustion (e.g., steam used to atomize fuel); pollutant control 
processes resulting in the net in-flow of water; and/or soot blowing processes usmg water to clean 
heat exchanger surfaces (commonly used in coal-fired systems). These sources of working fluid 
are in addition to: water formed from the combustion of hydrocarbon fuels; free water bom by the 



fiiel; and moisture carried by combustion air including air leakage. All such sources of effluent 
water are addressed by this invention through combustion stoichiometrics in combination with an 
ability to correct effluent data through use of optimization procedures. 

[0016] This invention adds to the technology associated with Input/Loss methods. 

Specifically The Input/Loss Method has been applied through computer software, installable on a 
personal computer termed a "Calculational Engine", and has been demonstrated as being highly 
usefiil to the operators of fossil-fired systems. The Calculational Engine receives data from the 
system's data acquisition devices. The Calculational Engine's software consists of the EX-FOSS, 
FUEL and HEATRATE programs described in '994 and '429, and in FIG.2 herein, and the ERR- 
CALC program described in '879 (also employed in '035 and '932) and in FIG.3 herein. ERR- 
CALC now incorporates the teachings of this invention. The Calculational Engine continuously 
monitors system efficiency on-line, i.e., in essentially real-time, as long as the thermal system is 
burning fiiel. The application of this invention to The Input/Loss Method significantly enhances the 
system operator's ability to predict tube failures and reduce outage time required for repair. 
[0017] The present invention provides a procedure for determining tube leaks in a recovery 

boiler using combustion stoichiometrics in combination with an ability to correct effluent data such 
that consistent fiiel chemistry is computed. 

[0018] The present invention teaches the mechanism of how a tube failure has been 

detected stoichiometrically, such detection being important to the system operator. Also, the 
present invention teaches how the location of a failed tube may be determined. 
[0019] Other objects and advantages of the present invention will become apparent when its 

general methods are considered in conjunction with the accompanying dravmigs and the related 
inventions of '994, '429, '563, '879,'035 and '932. 

[0020] The principles of this invention have been reduced to practice and installed for 

demonstration at two conventional power plants to determine the operability and fimctionality of 

this invention. These demonstrations have produced outstanding results demonstrating several 

identified tube failures. Furthermore, testing at a power plant burning PRB fiiel similar to black 

liquor in its fiiel water content, in which water soot blowing flow was removed fi'om the 

Calculational Engine's input (b^ = 0.0) which then computed an emulated tube leakage (bj > 0.0) 

with outstanding agreement to the measured soot blowing flow; see FIG.4. 

/// 

/// 

/// 

/// 

/// 



BRIEF DESCRIPTION OF THE DRAWINGS 



[0021] FIGURE 1 is a schematic representation of a recovery boiler illustrating the 

application of stoichiometric relationships, and also contains definitions of some of the important 
terms used herein. 

[0022] FIGURES 2 A and 2B is a block diagram of the general interactions and functions of 

the computer programs ERR-CALC, FUEL, EX-FOSS and HEATRATE used to implement this 
invention; herein collectively referred to as FIG.2. FIG.2 illustrates the "Fuel Iterations" involving 
FUEL, EX-FOSS and HEATRATE. 

[0023] FIGURE 3 is a block diagram of the principal functions of the error analysis 

computer program ERR-CALC which determines optimized Choice Operating Parameters, one of 
which is the tube failure flow rate incorporating the teachings herein. 

[0024] FIG.2 and FIG.3 relate the mteractions of the computer programs which implement 

this invention. The majority of the teachings of this invention are implemented in the ERR-CALC 
program. The FUEL, EX-FOSS and HEATRATE programs remain substantially unchanged as are 
used to implement the teachmgs of '994 and '429. However the FUEL, EX-FOSS and 
HEATRATE programs employ the results from ERR-CALC, including its calculated tube leakage 
flow rate, and thus assess the impact of such leakage on the thermal system in terms of boiler 
efficiency, fuel flow and system efficiency. Further, through energy balances on the steam 
generator's working fluid, and use of iterative procedures involving these programs, determination 
is then made within HEATRATE as to which heat exchanger within the steam generator contains 
the failed tube. 

[0025] FIGURE 4 is a plot of a sensitivity test of this invention. This test employed water 

soot blowing flow, entering the combustion space of a fossil-fired steam generator, to emulate a 
tube failure. The soot blowing flow, normally an input to the Calculational Engine, was left an 
unknown to be computed by the methods of this invention as a tube failure. 

DESCRIPTION OF THE PREFERRED EMBODIMENT 

[0026] To assure an appropriate teaching of this invention, its description is divided by sub- 

sections. The first two present nomenclature, definitions of equation terms, typical units of 
measure, and meaning of terms used herein (such as Choice Operating Parameters and System 
Effect Parameters), encompassing Paragraphs 0027 through 0043. The remaining sub-sections, 
representing the bulk of the teachings, are divided into four general groups: 
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1) the first group presents system stoichiometrics applied to recovery boilers and the 
determination of fuel chemistry based on effluents, these teachings support all subsequent 
disclosures herein (encompassing Paragraphs 0044 through 0062, employing equations 
numbered less than one-hundred); 

2) the next group presents the determination of boiler eflSciency for a black liquor-fired 
boiler as it influences both tube failure flow rate and determining tube failure location 
(encompassing Paragraphs 0063 through 0075, employing equations in the one-hundreds); 

3) the next group teaches how a tube failure may be detected based on an ability to correct 
Choice Operating Parameters using multidimensional minimization techniques, this ability 
being dependent, in part, on system stoichiometrics, the computed fiiel chemistry and boiler 
efficiency (encompassing Paragraphs 0076 through 0108, employing equations numbered 
in the two-hundreds); and 

4) the last group teaches how both the tube leakage flow rate and its location in the steam 
generator are determined using, as a foundation, the preceding teachings (encompassing 
Paragraphs 0109 through 0119, employing equations numbered in the three-hundreds). 

The remaining paragraphs present a conclusion, THE DRAWINGS and related teachings. 
Teachings of multidimensional minimization techniques, as directly applicable to this invention are 
also presented in '879, '035 and '932. The present invention expands the utility of Input/Loss 
methods to recovery boilers, and specifically builds upon and expands the utility of The Input/Loss 
Method described herein and in '994, '429, '879, '035 and '932, and in '563 as it teaches the L 
Factor. The methods described in '563 teach the foundations of the L Factor used in 
multidimensional minimization techniques; the L Factor is fiirther expanded as taught herein to 
encompass black liquor fuels used in recovery boilers. 

Definitions of Equation Terms with Tvpical Units of Measure 
[0027] Stoichiometric Terms: 

a = Moles of combustion O2 input to the system; moles/base. 
aP = O2 entering with system air leakage (typically via the air pre-heater); mole/base. 
^DRY-theor ^ Molcs of combustion O2 input to the system required for theoretical 
combustion associated with Dry (water fi*ee) fuel; moles/base. 
^Act ^ Concentration of O2 in combustion air local to (and entering) the system; molar flection. 
= Moisture in the entering combustion air; directly proportional to the ambient 
air's specific humidity, o) (Ib-water/lb-dry air); moles^ase. 

= (oa (1 .0 + <|)Act) Ndry-AIR / Nh20 
hpfi = Moisture entering with system air leakage; moles/base. 
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bj = Moles of tube leakage; i.e., water in-leakage entering and mixing with 
the combustion gases from leaks in heat exchangers; moles/base. 

b^ = Moles of known water in-leakage entering and mixing with the combustion 
gases not related to heat exchanger leaks (bj); moles/base. 
^Act ^ Total effluent CO2 at the system's boundary (i.e., Stack); moieis/base. 

g = Calculational effluent O2 at the system's boundary associated with a hypothetically 

segregated organic component of the black liquor, without air leakage; moles/base. 
^Act ^ Total effluent oxygen at the system's boundary burning black liquor fuel; moles/base. 
= (g-h 20^X^8 + aNkBL/2) + aP 
j = Calculational effluent H2O at the system's boundary, without air leakage; moles^ase. 
^Act ^ Total effluent water at the system's boundary (j + b^P); moles/base, 
^theor ^ Total cfflucnt water at the boundary based on theoretical combustion; moles/base. 
kgL = Effluent SO2 at the system's boundary from the reaction of Na2S04 
with CO2 obtained at the system boundary; moles/base, 
n^ = Molar quantities of dry gaseous effluents of combustion at the system boundary without 
air leakage; specifically those products associated with the following quantities: 
Daci^ (g-^2aNTNs+aNkBL/2), h, e^^, m^^, o^kBL and (1.0-aN)xa6; 
note: ^n^ =100 moles of dry gaseous effluent at the Stack is the assumed 
calculational "base" for Eq.(19BL), see FIG.l; moles/base, 
nji = Molar quantities of non-gas products of combustion at the system boundary without 

moisture associated with air leakage, specifically those products associated with the 
following quantities: j, xa^o, o^(T^o + ket), o^T^^q, c^nTna^ ^nTns» <^nTno v; 
see FIG.l and Eq.(19BL); moles^ase. 
Nj^ = Molecular weight of compound k. 

Rp^^ = Ratio of moles of dry gas from the combustion process before entering the air 
pre-heater to the diluted gas leaving, typically: (Moles of CO2 entering 
the air pre-heater) / (Moles of CO2 leaving the air pre-heater), defined 
as the Air Pre-Heater Leakage Factor; molar fraction. 
= (Moles of O2 entering the air pre-heater) / (Moles of O2 leaving the air pre-heater). 

Tro ^ Moles of effluent K2CO3 if burning black liquor (o^ = 1.0); moles/base. 

T|^^ = Moles of effluent Na2S04 if burning black liquor (% = LO); moles/base. 

T^Q = Moles of effluent NaCl if burning black liquor (o^ = 1.0); moles/base. 

T|^s = Moles of effluent Na2S if burning black liquor (o^ = 1 .0); moles/base. 
(Tno ^bl) ^ Moles of effluent Na2C03 if burning black liquor (a^ = 1.0); moles/base. 



V = Moles of effluent free carbon, or its equivalence, typically found 

in the recovery boiler's smelt; moles/base. 
X = Moles of As-fired fuel required per 100 moles of diy gaseous effluent; moles/base, 
^eor ^ Moles of As-Fired fuel associated with theoretical combustion; moles/base. 
^DRY-theor ^ Molcs of Dry fiicl associatcd with theoretical combustion; moles/base. 
^MAF-theor ~ Molcs of Moisture-Ash-Frec fiiel associated with theoretical 
combustion; moles^ase. 
xttjQ = Inert matter in As-Fired fiiel, the terms "inert" and "ash" 
are used interchangeably; mole-inertl^ase. 
a^r = As-Fired (wet-base) fuel chemistry constituent k per mole of fuel; 

Y,a^ = LO, where: k= 1,2,3,4,5,6,10,14,15,16; see Eq.(19BL) 
therein for terms; mole-k/mole-fuel. 
^MAF-k ^ Moisture- Ash-Free (MAF) fiiel constituent k per mole of MAF fuel; 

Z^MAF-k ^ 1-0' where: k = 1,3,4,5,6,14,15,16; see Eq.(19BL) 
therein for terms; mole-k/mole-fuel. 
P = Air Pre-Heater Dilution Factor (ratio of air leakage to true combustion air); molar fraction 

P MOO(^A^ -L0)/[ai?Act(1.0 + (t)Act)] 

= Kronecker function: unity if black liquor (bearing at least sodium) is being 

employed in the fuel, otherwise zero; unitless. 

^p^^^ = Ratio of non-oxygen gases (Nj and Ar) to oxygen in the combustion air; molar ratio. 

c|)Act^(1.0- Act)/^Act 

<j)Rgf = Reference ratio of non-oxygen gases (principally N2 and Ar) to oxygen in the 

combustion air, taken as 3.7737245; molar ratio. 
<Pr = Reduction Efficiency; see Eq,(19BL); molar ratio. 

- / (Tna Tns) 

(ps = Sulfur to Smelt ratio; see Eq.(19BL); molar ratio. 

- (Tna + Tns) / [ONXttio + On(Tno + kBL) + o^T^o + Tna + + Tnc) + v] 
<Pn = Sodium to Carbon ratio in the fiiel; see Eq.(19BL); molar ratio. 

- ^MAF-14 / ^MAF^ 

[0028] Multidimensional Minimization Terms: 

F (x) = Objective function, a fimctional relationship using the independent variables 3c ; unitless. 
/() => Indicates a general functional relationship; for example, the expression: 

HHVj^ = /[fuel chemistry(A)], means that HHV^ is a fiinction of fiiel 
chemistry (which in-tum is a function of the vector A). 
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C i = Correction factor to be applied to an initial Choice Operating Parameter i; 
see Eqs.(21 IS) through (220) for nomenclature, e.g., Cj^ is the 
correction factor for Stack CO2 as referring to Eq.(21 1 S), etc.; unitless. 
HHVj^ = Higher heating value as used by the minimization techniques as a System Effect 
Parameter, here subscript k3 refers to either an uncorrected As-Fired (HHV^), 
Dry or MAP heating value; Btu/lbm^, Btu/lbrnQRY Btu/lbmj^. 
HHV|^.Rgf = Higher heating value used as a Reference System Effect Parameter; 
Btu/lbm^, Btu/lbm^RY Btu/lbmj^. 
Jq = Bessel function of the first kind of order zero. 
= Bessel function of the first kind of order one. 
Lqq2 = L Factor for dry effluent CO2 used by the minimization techniques 

as a System Effect Parameter; (Ibm effluent C02)/million-BtUpye,. 
^Vuei ^ L Factor for dry gaseous effluents used by the minimization techniques 

as a System Effect Parameter; (Ibm dry gaseous effluent)/million-BtUpugj. 
L^i = Generic L Factor, here subscript kl refers to either L'p^g, or Z'C02- 
^C02-Ref ^ L Factor for dry effluent CO2 used as a Reference System Effect 

Parameter; (Ibm effluent C02)/million-BtUpyei. 
^Vuei-Ref ^ L Factor for dry gaseous effluents used as a Reference System Effect 
Parameter; (Ibm dry gaseous effluent)/million-BtUpyg,. 
myy; = Fucl flow ratc, an As-Fired mass flow quantity (i.e., wet with water and fuel 

minend matter), as may be computed by Input/Loss methods; also may be used 
by minimization techniques as a System Effect Parameter; Ibm^/hour. 
^AF-PLT ^ Th^ system's measured fuel flow, an As-Fired quantity (i.e., wet with water and fiiel 
mineral matter), also termed the system's "indicated fuel flow"; also may 
be used as a Reference System Effect Parameter; Ibm^/hour. 
m-p = Tube leakage flow rate; i.e., mass flow rate of water in-leakage entering the combustion 
gas path from leaks in a heat exchanger, a Choice Operating Parameter; Ibm/hour, 
Ml = Dilution Factor applied to System Effect Parameter L^^\ Ml > 0.0; unitless. 
M^ = Dilution Factor applied to System Effect Parameter m^; M^ > 0.0; unifless. 
Mji = Dilution Factor applied to System Effect Parameter HHV,^; M^ >0.0; unitless. 
My = Dilution Factor applied to System Effect Parameter WFy^Q, My > 0.0; unitless. 
Sj = Scaling factor for the independent variable Xj ; reciprocal units of measure of Aj . 
Sj = Pre-scaling factor used to adjust S|; unitless. 
WF^ = As-Fired fiiel mass fraction of substance j; fraction. 
WFyqq= As-Fired fiiel water mass fraction (also termed WF^, used as 
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a System Effect Parameter; fraction. 
^^H20-Ref Reference As-Fired fuel water mass fraction used by the 

minimization techniques as a Reference System Effect Parameter; fraction, 
jc = Vector of independent variables, x = (x^, x^, ... ), as based on scaled Choice Operating 
Parameters (not to be confused with the term for moles of As-fired fiiel, x); unitless. 
Aj = Choice Operating Parameter i, see the specific parameter for units of measure, 

and Eqs.(21 IS) through (220) for defmitions. 
A = Vector of Choice Operating Parameters, which is user selected; for example, one 

selection might include: A = (A^^, A2S, A3, A^, Ayg, A^); seeEqs.(211S)thru(220). 
Aq.! = Initial Choice Operating Parameter i, before application of a minimization technique 
and based on the system's raw instrumentation signal, a previous converged 
solution, an estimate, and/or as otherwise obtained. 
Ap.i = Converged (fmal) Choice Operating Parameter i, after application of a niiiiiniization 

technique to Aq.j, and, thus corrected and applicable to all system thermal analyses. 

[0029] Quantities Related to System Terms: 

AF = Air/Fuel ratio defined by indicated air flow and m^y^.pLj; unitless mass ratio. 
BBTC = Energy flow to the working fluid from combustion gases (w/o tube leakage); Btu/hr. 

Cp = Heat capacity; Btu/lbm-R. 

h j = Specific enthalpy of substance j; Btu/lbm. 
HBC = Firing Correction; Btu/lbmyy7. 

HHVP = As-Fired higher heating value, based on HHV^ and used in system 

evaluations as corrected for a constant pressure process; Btu/lbm^y?. 
HNSL - Non-Chemistry & Sensible Heat Losses; Btu/lbm^. 

HPR^j^ = Enthalpy of Products from actual combustion (HHV- or LHV-based); Btu/lbm^. 
HPRideai " Enthalpy of Products from ideal combustion associated with an 

obtained heating value at T^^ (HHV- or LHV-based); Btu/lbm^^. 
H^Act = Enthalpy of Reactants associated with actual firing conditions 

(HHV- or LHV-based); Btu/lbm^. 
HRX^ai = Enthalpy of Reactants associated with an obtained heating value 

at Tcaj (HHV- or LHV-based); Btu/lbm^. 

HR = System heat rate (HHV-based,/«HHv; or LHV-based, /«LHv); Btu/kWh. 
HSL = Stack Losses (HHV- or LHV-based); Btu/lbm^. 
LHVyyr = Lower heating value based on the measured or calculated 
higher heating value (HHV^^); Btu/lbm^^:. 
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LHVP = As-Fired lower heating value, based on LHVy^ and used in system 

evaluations as corrected for a constant pressure process; Btu/lbm^y?. 
m-pAh =Energy flow from tube leakage; Ah is the enthalpy difference between the last 
heat exchanger effected by the leakage (typically the Final Superheater or 
Reheater) and the heat exchanger producing the leakage; Btu/hr. 
= Ambient pressure local to the system, psiA. 
QsAH ^ Energy flow delivered to steam/air heaters, Btu/hr. 
T = Temperature; F. 

= Ambient temperature local to the system, F. 
Tcai - Calorimetric temperature to which heating value is referenced, F. 
Tstack " Boundary temperature of the system effluents, defines the "Stack"; F. 

= Brake power associated with in-flow stream fans (e.g.. Forced Draft fans); Btu/hr. 
= Brake power associated with out-flow stream fans (e.g., Induced Draft fans); Btu/hr. 
^output ^ G^'oss power generated from a power plant; kWe. 
T]sYs ^ System efficiency (HHV-based, t|sys-hhv5 LHV-based, rj^Ys-LKv); unitless 
tiq Boiler efficiency (HHV-based, rig^j^^y; or LHV-based, t|3.lhv); unitless. 
r]Q = Combustion efficiency (HHV-based, Tj^^^H^y; or LHV-based, r|c.LHv); unitless. 
r\p^ = Boiler absorption efficiency; unitless. 

[0030] Subscripts and Abbreviations: 

Act = Actual value obtained from the operating thermal system. 

AF = As-Fired fuel at the thermodynamic boundary (i.e., wet with water and inert matter). 
BL = Black liquor fuel. 
DRY= Dry chemical base (i.e., free of water). 
MAF = Moisture- Ash-Free chemical base (i.e., free of water and free of inert matter). 
Ref = Reference value. 
T = Tube failure associated with a heat exchanger leakage, 
theor = Refers to conditions associated with theoretical combustion. 

Meaning of Terms 

[0031] The words "Operating Parameters", as taken within the general scope and spirit of 

the present invention, mean common data obtained from a thermal system applicable to the 
thermodynamic understanding of that system. The following quantities may be included in the 
definition of Operating Parameters, they are not encompassing but considered typical of a 
minimum set of data required for such thermodynamic understanding: effluent CO2, O2, and SO2 
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concentrations determined at the Stack, or before the air pre-heater (Boiler side of the air pre- 
heater); the mass, wet-base ratio of the indicated combustion air flow at the system's fuel 
combustors, to the system's indicated fiiel flow, termed AFj^^ (note that AFj^^^ is obtained only for 
the determination of inerts as taught in '994); effluent H2O concentration measurement, or 
assumptions made (or as otherwise may be determined); effluent temperature measurement, that is 
the average temperature associated with the combustion gases at the system boundary (caution 
must be exercised in measuring non-stratified gas flows); the inlet/outlet ratio of CO2 (producing 
R^^^ as is preferred), CO or O2 (producing /? Act) across the air pre-heater where these ratios could 
be obtained on-line, off-line, based on periodic testing or judgement which are used for the 
determination of air pre-heater leakage; determination of fuel temperature at an appropriate system 
boundary; air psychrometric measurements, or as otherwise determined, at the system boundary 
(e.g., dry and wet bulb temperatures, or dry bulb and relative humidity, or dry bulb and dew point 
temperatures); quantities comprising the system's Firing Correction term, HBC; the discharge 
temperatures of the air as it exits each air heating or cooling device but before it reacts with the fuel 
(for example, such devices might include the air pre-heater, forced-draft fan, steam-to-air heater, 
etc.); the total energy flow deposition to the working fluid from the combustion gases (in typical 
units of measure of Btu/hr); and similar quantities. Regards the total energy flow deposition, for a 
typical recovery boiler, such measurements typically include feedwater flow to the steam 
generator, feedwater pressure and temperature, determination of the steam flow from the steam 
generator if different than the feedwater flow, steam pressure, steam temperature or quality (or 
assumed quality), and, if applicable, reheat flows, and reheat inlet and outlet pressures and 
temperatures. If employing a Reheater heat exchanger, determination of accurate reheat flows 
generally requires understanding of steam turbine flow distributions (involving high pressure 
turbine shaft seals, steam flows to feedwater heaters, turbine bypass leakages, attemperation spray 
flows and the like). The total realizable energy flow (BBTC - m^Ah) includes the eflFects of tube 
leakage flow rate, determined by this invention, as assigned to a particular heat exchanger. 
Specifically, when having determined there is a tube leak, the leaking heat exchanger's outlet flow 
of working fluid is reduced by the tube leakage, thus reducing the energy flow BBTC by mjAh. 
[0032] The words "Choice Operating Parameters", as taken within the general scope and 

spirit of the present invention, mean any sub-set of Operating Parameters which directly impact 
system stoichiometrics, and thus may impact the determination of fuel chemistry. This invention 
assumes that Choice Operating Parameters may have error, said error may adversely effect the 
determination of fuel chemistry, but said error may be corrected as taught herein through 
optimization methods. In the Preferred Embodiment Choice Operating Parameters are selected by 
the user of this invention from an available set This available set of Choice Operating Parameters 
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includes the following ten: 1) effluent CO2 concentration measured at the Stack or Boiler; 2) H2O 
concentration measured, or as otherwise may be determined, at the Stack or Boiler; 3) the mass, 
wet-base ratio of the indicated combustion air flow at the system's fuel combustors, to the system's 
indicated fiiel flow, the Air/Fuel ratio termed AF^^{, 4) the Air Pre-Heater Leakage Factor, termed 
Rj^^^; 5) the concentration of O2 in the combustion air local to the system, or as otherwise 
determined, termed Aj^^ (leading to the determination of (|)Act); 6) the system's indicated limestone 
mass flow rate, termed m^^; 7) effluent O2 concentration measured at the Stack or Boiler; 8) mass 
flow rate associated with a heat exchanger tube leakage flow rate, termed m^; 9) Sodium to Carbon 
ratio found in the fuel ((p^); and 10) the Reduction Efficiency ((p^). 

[0033] The words "Reference Fuel Characteristics", as taken within the general scope and 

spirit of the present invention, mean the identification and development of intrinsic chemical 
relationships which allow resolution of fuel chemistry. Discussions of Reference Fuel 
Characteristics are presented in Paragraphs 0048 through 0058. For the Preferred Embodiment, 
Reference Fuel Characteristics include an average or typical fiiel chemistry and associated MAF 
heating value, preferably based on historical data collections of ultimate analyses of the fuel's 
elemental composition, with fiiel water and fuel inerts (typically reported as weight fi'actions 
summing to unity, leading to % molar fi'actions), herein termed fiiel chemistry constituents. If 
ultimate analyses are not available they may be approximated fi'om proximate analyses as taught in 
Gill's text (cited below, not preferred). Reference Fuel Characteristics typically include a fuel 
hydrogen versus fiiel carbon relationship, typically a MAF relationship; a relationship based on 
historical data Reference Fuel Characteristics may also include relationships of: MAF oxygen 
versus MAF carbon; MAF sulfur versus MAF carbon; MAF potassium versus MAF sodium; and 
MAF chlorine versus MAF sodium. Reference Fuel Characteristics may also include the computed 
values of Z pug|.j^gf and ^co2-Ref» computed using the reference fiiel chemistry. For black liquor fiiel 
Reference Fuel Characteristics may also include the variables (p^, (pg and (p^- Reference Fuel 
Characteristics may also include whether the variability of fiiel water and fuel inert flections in the 
As-Fired condition is predictable, or not. For any given fuel: fuel water may be held constant 
(including zero); fuel inerts may be held constant (including zero); fimctionalities may be observed 
and applied, for example, aMAF-10 ~/ (HH^mafX and/or fuel water and/or fuel inerts may be 
treated as unknowns). The Preferred Embodiment for black liquor fiiels is to treat fuel water as an 
unknown, and fiiel inerts as a constant. All of these possible variations for the treatment of fiiel 
carbon, fuel sodium, water and inerts may be included as a portion of Reference Fuel 
Characteristics. Reference Fuel Characteristics may also contain fitting constants associated with 
all correlations relating dependent fuel quantities to independent fiiel quantities, typically the 
principle independent fuel quantity is fiiel carbon. When applying the Preferred Embodiment of 
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this invention, Reference Fuel Characteristics also contain reasonability limits (i.e., numerical 
minimum and maximum limits) of the following: computed fuel elemental constituents; fuel water 
fraction; fiiel inert fraction; the correction factors being applied to initial Choice Operating 
Parameters; the maximum rates of change of these parameters; and the maximum rates of change 
of any or all of the selected Choice Operating Parameters. Such minimum and maximum limits are 
preferably based on engineering judgement supported by historical data collections of ultimate 
analyses of the fiiel, historical experience of instrumentation on which the selected Choice 
Operating Parameters are based, and historical records of computed correction factors to the initial 
Choice Operating Parameters, termed Aq.j. 

[0034] The words "System Effect Parameters", as taken within the general scope and spirit 

of the present invention, mean any parameter of the thermal system or its fuel which dh-ectly 
impact the determination of system efficiency. In the most general sense System Effect Parameters 
include any parameter used in Eq.(331A), (33 IB), (332A) or (332B) which compute system heat 
rate and thus system eflBciency. For the Preferred Embodiment, System Effect Parameters include 
the following four types of quantities: the L Factor (either I p^gj or Lqqj)\ the computed As-Fired 
fuel flow rate (m^); the higher heating value (either HHV^, HHVj)ry or HHVmaf); and the As- 
Fh-ed fiiel water faction {WFy^^ which may be used to convert HHVdry to HHV^. The 
computed L Factor effects fuel chemistry which effects heatmg value and boiler efficiency, thus 
has an unmediate impact on system efficiency. "Reference System Effect Parameters" are constant 
and targeted (i.e., desired) System Effect Parameters to which the System Effect Parameters are 
numerically driven by the minimization techniques through optimizing a selection of Choice 
Operating Parameters. 

[0035] The words "Input/Loss methods", as taken within the general scope and spirit of the 

present invention, mean any method or combination of methods in which one or more of the 
following parameters is determined based on Operating Parameters and a selection of Choice 
Operating Parameters: fuel flow, effluent flow, emission rates, ftiel chemistry, fuel heating value, 
boiler efficiency, and/or system efficiency. In addition to '994, '429, '563, '879, '035 and '932, 
and their related provisional patent applications and Continuation-In-Parts, Input/Loss methods 
include the methods of '470 and '420. The words "The Input/Loss Method" refers specifically to 
the collection of technologies described in '994, '429, '879, '035 and '932, and in '563 as it teaches 
the L Factor, and to any relevant provisional patent applications and Continuation-In-Parts, in 
addition to the teachings disclosed herein. 

[0036] As used herein, the words "Calculational Engine" refers to a computer in which 

software descriptive of The Input/Loss Method is installed. 
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[0037] As used herein, if used, the words "obtain", "obtained", "obtaining", "determine", 

"determined", "determining", "determination", "establish", "established" or "establishing" are 
defined as measuring, calculating, computing, assuming, estimating or gathering from a database. 
[0038] As used herein, the words "monitoring" or "monitored" are meant to encompass 

both on-line monitoring (i.e., processing system data in real time) and off-line monitoring (i.e., 
computations involving static data). A "monitoring cycle" is meant one execution of the processes 
described in FIG. 2, which encompasses FIG. 3. 

[0039] As used herein, the meaning of the words "smoke Stack" or "Stack" or "system 

boundary" are defined as the physical boundary of the thermal system where gaseous combustion 
effluents exit, entering the local environment; refer to 42 in FIG.l, further discussed within THE 
DRAWINGS. Solid effluents not leaving the Stack (e.g., molten smelt 37) are referenced to the 
generic system's boundary 44 in FIG.l. 

[0040] As used herein, the meaning of the words "Boiler" or "Boiler Effluent" are defined 

as the region 35 in FIG.l, or generically between the physical exit of the system's region 34 in 

FIG. 1 and entrance to its air pre-heater 36 in FIG. 1 ; see THE DRAWINGS. 

[0041] As used herein, the meaning of the words "Fuel Iterations", are defined in 

conjunction with a detailed description of FIG.2 found within THE DRAWINGS. 

[0042] As used herein, the meaning of the word "indicated" when used in the context of 

data originating fi-om the thermal system is defmed as the system's actual and uncorrected 

measurements of a physical process (e.g., pressure, temperature, mass flow, volumetric flow, 

density, and the like) whose accuracy or inaccuracy is not assumed. As examples, a system's 

"indicated fuel flow" or its "indicated limestone flow" denote system measurements the accuracy 

of which is unknown (they are "as-is", with no judgement applied). Such indicated measurements 

are said to be either correctable or not. If not correctable, it may be that the associated computed 

value from Input/Loss methods tracks the indicated value over time (the indicated not being 

corrected per se). In the case of indicated limestone flow when used as a Choice Operating 

Parameter (A^), it is directly corrected as taught by this invention. In the case of indicated fuel 

flow when used as a System Effect Parameter, it may be shown that the computed fiiel flow, m^^?, 

tracks the indicated fiiel flow, m^^.pLT, through adjustment of the Dilution Factor M^, 

[0043] As used herein, the meaning of the words "Reduction Efficiency" is defmed by the 

variable (p^; the meaning of "Sulfur to Smelt ratio" is defined by the variable cp^; and the meaning 

of "Sodium to Carbon ratio" in the fuel by the variable (p^. 

/// 

/// 

/// 
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System Stoichiometrics 

[0044] Any study of the combustion of fossil fuels necessitates the formulation and use of a 

combustion equation. Combustion equations used by several Input/Loss methods are described in 
^470, '420 and '994 by their designated Eq.(29), in '429 by its Eq.(19), in '879 by its Eq.(19-corrX 
and in '035 and '932 by their Eq.(19B). This invention's methods are taught through a combustion 
equation defined by Eq.(19BL) herein. Through Eq.(19BL) stoichiometric terms are self-defined, 
including the use of the hj term important to this invention and inclusion of fuel components 
typically found in recovery boilers (hydrocarbons, sodium, potassium and chlorine). Eq.(19BL)'s 
nomenclature is unique in that brackets are used for clarity: for example, the expression 
"xa2[H20]" means the moles of fiiel water/base, algebraically simply xa2; the expression 
"o^xai^\Na2T means the moles of foel diatomic sodium/base, algebraically simply ONxai4; the 
expression "(d - o^k^i)j^^[C02Y' means the effluent moles of carbon dioxide/base, algebraically 
simply (d - ONk3L)Act which defines the term D^^.^ ; etc. The stoichiometric base of Eq.(19BL) is 
100 moles of dry Stack gas (i.e., at the thermodynamic boundary). Also, note that in combustion 
equations such as Eq.(19BL) describing commercial recovery boilers, that the assumption is 
typically made of possible extraneous water leakage into, and mixing with, the products of 
combustion. The flow of such water is assumed known (either measured or otherwise reasonably 
estimated) and is denoted by the symbol b^. It is common industrial practice in the United States to 
use this symbol, b^, to denote the quantity of known water in-leakage such as steam used to 
atomize fuel, water used for soot blowing, etc. Such b^ in-leakage is apart from water formed from 
combustion of hydrocarbon fuels, tiie xa5 term; apart from free water bom by the fiiel, xa2; apart 
from moisture carried by the combustion air and air pre-heater leakage, (1.0 + and apart from 
tube leakage, the b^ term. 

X [ ai[N2] + a2[H20] + a3[02] + a4[C] + a^\H2\ + ^^[S] + ajoEInert] 
+ aNai4[Na2] + o^a^^\K2\ + aNai6[Cl] ] 

As-Fired Fuel 
^z\^2^]\n-Leakage WP20]Tube-Leakage 

+ [ (1.0 + P)(a[02] + a4)Aet[N2] + bAPjO]) U 
= (d - o^kBiJLCOj] + (g + 2onTns + ONkBL/2)[02] + h[N2] + j[H20] 

+ xaioPnert] + e^^JCO] + m^JNO] + o^T^q + k^Ol^a^CO^] 
+ OnTko[K2C03] + ONTNA[Na2S04] + ONTNsLNajS] 
+ aNTNc[NaCl] + ONkBL[S02] + (1.0 - aN)xa6[S02] 

+ v[CRefiJ + [ P(a[02] + a(j)Act[N2] + h^lHiO]) lAirLeakage (19BL) 
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[0045] Resolution of Eq.(19BL) is had when all n-^ and n^i quantities have been determined. 

Minor component terms of Eq.(19BL) are typically resolved either through direct measurement 
(e.g., for CO and NO), or assume zero values, or through obtained relationships. All Minor 
Components typically have only low parts-per-million concentrations thus have little impact An 
exception to this is the v term describing refuse carbon found in the smelt. The true importance 
and functionality of Eq.(19BL) to The Input/Loss Method lies in the fact that consistency of molar 
balances is needed for successful system understanding, for conservation of mass flows and for 
resolution of fiiel chemistry, these needed for detection of tube failures. For clarity the following 
major terms are associated with system stoichiometrics: 
Total effluent (boundary) water = 1^^,^ = j + b^p 
Boiler oxygen before air leakage (g Act) - ^Act(g ^o^Ti^g + a|^kBL/2) 
Total effluent (boundary) oxygen = G^ct (g 2a|^T|^s + 0|^kQL/2) + aP 
Total effluent (boundary) carbon dioxide = = (d - o^^kgL) 
Total effluents referenced to the boundary = Enj + Enjj + P(a + ac^^ct ^ b^) 
Total effluents before air leakage, referenced upstream of the air pre-heater 

Dry combustion air without air leakage referenced to the boundary = (a + a<l)Act) 
Wet combustion air without air leakage referenced to the boundary = (a + ^^^ct '^a) 
Dry air from air leakage found at the boundary = p(a + ^^Act) 
Total wet combustion air and air leakage found at the boundary 
= (1.0 + P)(a-fa(l)Act + bA). 

[0046] Eq.(19BL) is unique in describing at least three features of critical importance when 

determining fuel chemistry using one of the Input/Loss methods. The critical features include: 1) 
its ability to address air pre-heater leakage through application of the Air Pre-Heater Leakage 
Factor, /?Act' through the Air Pre-Heater Dilution Factor, p; 2) the ability to describe effluent 
concentrations on either side of the air pre-heater, again through application of i?Act' 3) the use 
of an explicit ^p^^ term allovsdng for variable O2 concentration in the system's local combustion air. 
Air pre-heater leakage dilutes all combustion effluents with moist air from the local environment, 
thus all important effluents H2O, CO2 and O2 used for this invention are altered. Furthermore, 
many times, although not always, a power plant's more precise effluent measurements may be 
found on the air pre-heater's inlet (economizer outlet or Boiler), and not at the air heater outlet (or 
Stack); refer to FIG.l. Although most environmental regulation requires effluent measurements at 
the system's boundary, translation between the air heater inlet and outlet measurements is many 
times essential. Eq.(19BL) allows for such translation through the i?Act term, defined above such 
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that 100 moles of dry gas are computed both at the upstream and downstream locations of the air 
pre-heater; see "Boiler" of FIG.l. Thus effluents may be used by the present invention either 
upstream or downstream of the air pre-heater; refer to the G^^t ^Act terms defined above, 
allowing conversion between measurements with and without air leakage. For example, 
combustion gas conditions for oxygen and water upstream of the air pre-heater and after exiting the 
heat exchangers and combustion region, see FIG.l, would employ the terms: 
(g+2oj^TNs-i-0|^kBL/2)i?y^ct j^Acf That is, one would actually measure 

(g+20|^T|sjs+a|sjkBL/2)/?^^^ moles of dry O2 upstream of the air pre-heater and after exiting the heat 
exchangers and combustion region as based on 100 moles of dry gas found at that location. As a 
fijrther example, the moles of dry air found upstream of the air pre-heater is given by: Rj^^^ a(1.0 + 
4>Act)- Combustion gases downstream of the air pre-heater typically exit the system to the 
environment (i.e.. Stack), in other words the gaseous effluent boundary of the system (100 moles of 
dry gas at the Stack includes air leakage). If limestone is injected into the combustion process to 
control effluent SO2 it will create additional effluent CO2; further, it could decrease the effluent 
H2O if the sulfate product is matrixed with water, CaS04-zH20. Thus such effects must be 
considered, and are taught in '994 referenced herein. Although the specific use of limestone is 
taught in '994 it is not included herein for clarity of presentation. Of course CO2 and H2O are two 
important effluents to the present invention. In addition to the basic stoichiometrics afforded, 
Eq.(19BL) allows numerous and obvious determinations of molar and mass ratios. 
[0047] Based on these teachings and those presented in '994, the following further explains 

the importance of the Air Pre-Heater Leakage Factor, Rj^^ and the Air Pre-Heater Dilution Factor, 
P, their defmitions and developments. Consider that air in-leakage associated with a fossil-fired 
system, and as commonly associated with in-leakage at the system's air pre-heater, is defined by 
the American Society of Mechanical Engineers' Performance Test Code 4.3 (1974) as the moist air 
leakage divided by the wet combustion gas. Typically the wet combustion gas is taken at the gas 
inlet of the air pre-heater (i.e.. Boiler, or economizer outlet before the air pre-heater). That is, as 
defined herein using Eq.(19BL) nomenclature, noting that 100 moles of dry gas is the bases at the 
Boiler, is given by: 

Wet APH Leakage = ^Act-£(9iiiAct-i^Ai^MoistAir- (20) 

(100 + i?Actj)NwetGas 

where, as defined above: 

R^^ = (Moles of CO2 entering the air pre-heater)/(Moles of CO2 leaving the air pre-heater). 
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The expression for /J^ct is equivalent to [Moles of Boiler COj] divided by [Moles of Stack COj]. 
The Air Pre-Heater Dilution Factor is then developed by performing a total dry gaseous effluent 
molar balance at the Stack: 

100 moles dry gaseous effluent at Stack = Enj + P(a + a<|)Act) (21) 

then solving for P: p = (100 - En{) / (& + a^^^^). The stoichiometric base of Eq.(19BL) implies 
that 100 moles of dry gaseous effluent upstream of the air pre-heater (Boiler) is given by Rj^c^ni, 
see FIG. 1; therefore: 

P=(100 - 100//?Act)/[a(1.0 + (|)Act)] 
- 100(/?A^ - 1.0) / [/?Act a(1.0 + <|)Aet)]. (22) 

If, instead of obtaining the ratio of COj across the air pre-heater, the ratio of O2 is obtained and, 
following the teachings herein and those of '994, the following may then be developed: 

= (Moles of O2 entering the air pre-heater) / (Moles of leaving the air pre-heater). 
where, converting from R'^ to Rj^^ using algebraic manipulations results in, when measuring 
Stack O2: 

ji^^ , 100 -ig'ActGAct(1.0 + <|>A^ .23) 

100-GAet(1.0 + <|)Act) 

If measuring Boiler O2 (for Eq.(24) termed g\^: 

, -iOQgAct-l^ActglActf IrO + (|)ActX- (24) 
100^Act-gAct(1.0 + <|)Ac.) 

There are, of course, a number of variations to these formulations, such as employing 100 moles of 
wet effluents at the Stack, thus replacing Eq.(21) with: 

100 moles wet effluent at Stack = (En^ + j) + p(a + sa^j^^ + ) (25) 

or using an oxygen base for the wet effluents at the Stack, thus: (En, + J^cd^^ + 4)Act)^ 

using a combustion equation which is based on a mole of fiiel carbon (xa4); etc. What is important 
to this invention, important to The Input/Loss Method, and important to any of the Input/Loss 
methods, is that the Air Pre-Heater Leakage Factor (/?Act) allows gaseous measurements to be 
employed on either side of the system air in-leakage. Typically, but not always, is measured in 



-23- 

the combustion gas stream inlet to the air pre-heater (Boiler), while CO2 is measured at the Stack 
(downstream from the air pre-heater). 

[0048] Operators of recovery boilers will commonly determine certain intrinsic fuel and 

effluent ratios which assist in understanding the performance of the combustion process, i.e., 
performance ratios. Such ratios are determined through periodic fiiel analysis, and analysis of the 
combustion residuals (also termed "smelt" for recovery boilers). Smelt is defined as the molten 
inorganic compounds leaving the fiimace after combustion; using Eq.(19BL) nomenclature it 
consists of: xa^Q, ^^(Tisjo+kBL), o^Tj^O' ^n^na' ^n^ns^ ^n^nc Common performance 

ratios include: the Reduction Efficiency ((p^); the Sulfiir to Smelt ratio ((p^); and the Sodium to 
Carbon ratio in the fuel ((p^). Reduction Efficiency judges the ability of the combustion process to 
further reduce Na2S04. For theoretical combustion, (p^ = 0.0; for complete reduction of Na2S04, 
(Pr = 1.0. The Sulfur to Smelt ratio is employed in the Preferred Embodiment of this invention. 
The Sodium to Carbon ratio in the fuel is a key parameter, determined periodically by laboratory 
analysis, as it allows stoichiometric inter-dependency between the organic and inorganic fuel 
components. Other such parameters may be formed which are similar to the ratio. As taught 
herein, what is required is some defmed or otherwise obtained dependency between the organic and 
inorganic fuel compounds if fuel chemistry is to be determined based on gaseous effluent and smelt 
measurements. Although the Preferred Embodiment of this invention employs the (p^ ratio, on-line 
measurements of the As-Fu-ed fiiel and/or on-line measurements of the smelt could prove as useful 
as the (Pn ratio, and/or otherwise aid in its determination. 

[0049] After establishing system stoichiometrics, the next stage of the process involves the 

recognition that a given fuel has an unique chemical composition, thus when bumed will yield 
unique concentrations in its gaseous effluent. The gaseous effluent concentrations are used to 
compute fiiel chemistry, with this chemistry heating value and boiler efficiency are computed, in 
turn this information allows the detection of tube failures and their location. However, key to this 
invention is that its methods recognize that gaseous effluent concentrations are not accurate and 
require correction to achieve stoichiometric consistency. The gaseous effluents fi^om any fossil 
combustion process, including black liquor combustion, are N2, COj, H2O, O2 and SO2. HjO, 
when effluent fi-om combustion, is in its superheated phase, thus acting as a gas. The source of N2 
is principally the air used to bum the fuel and has little chemical reactiveness, thus its sensitivity to 
the fueFs chemical composition is not significant. SO2 effluent concentrations are generally in the 
parts per million thus its impact has minor importance. SO2 effluent is formed fi-om either Na2S04 
reduction or, if burning a pure hydrocarbon, directly from fiiel sulfur (1.0 - o^)xa^. 
[0050] As an intrinsic chemical relationship, the relative concentrations of carbon, a^, and 

hydrogen, a^, found in any fossil fuel will have significant impact on the relative concentrations of 
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CO2 and H2O found in the effluent. In addition, these effluents will be influenced by the 
following: O2 used to bum the fuel (i.e., the Air/Fuel ratio); fuel water, in-leakage of water 
including tube leaks; and water in the combustion air. This implies that the molar fractions of CO2, 
H2O and O2 present in the effluent (the system's boundary, i.e., its smoke stack or translated from 
air pre-heater inlet data) must be unique relative to the supplied fuel and supplied combustion air. 
[0051] The following elemental molar balances may be derived from the combustion 

equation, Eq.(19BL). The Fj^ expressions are simply convenient groupings of quantities, 
principally comprising measured effluents (known values) which have the greatest influence on the 
individual elements of interest. Minor fuel terms, aMAp-j? carried within expressions are 
multiplied, initially, by an estimated fuel moles, x^/^. These minor terms are quickly resolved 
when converging on x^^. Given these groupings, the Ty. expressions of Eqs.(36) through (41), 
with solution of the moles of combustion oxygen (the term "a") as discussed below, may be treated 
as known quantities. The elemental wet fuel components typically associated with black liquor 
fuels are considered unknowns, as are the fuel moles, these include the following: aj, a2, a^, a^, a^, 
ttg, ttjo, ^14, a|5, and x. Many recovery boilers use supplementary firing with gaseous fuel such 
as methane or fiiel oil. Inclusion of such fuels is taught in '994 being described by the combustion 
equation term: xaQ[CYRH2R] of that patent. Note that The Input/Loss Method accounts for all non- 
black liquor fuels, such as methane or fiiel oil, but they are assumed to have known chemistries 
with known fuel flows. Mixed fuels are input to the system simulator (e.g., EX-FOSS) as a 
composite, combining them by using the FUEL program, although the resultant computed fuel 
chemistry, after proper weighting, is assigned to the black liquor by assuming the other fuels have 
fixed or predictable, and known, chemistries. Alternatively, such assignment may be made to a 
combined subset of fuels. 



(30) 
(31) 



x(a5 -f a2)/2 = 



x(a3 + a2/2) = ro2 + 3(xMAFrNA + 



^MApaMAF^ (^N - 2aN<pR + 1 .0) 
xa4 = rco2 x^afFna + o^kpL 

xa,4 = aN(TNo + kBL) + ^N(TNA + TNs + TNc/2) 



(32) 
(33) 
(34) 
(35) 



where: Fn2 = 100- [DAct + GAet + eAet + niA,/2 

+ ^N^BL + (1.0 - aN)XMAFaMAF-6] 

- 1004)Act(^Act- l-0)/[i?Act(*Act+ 1-0)] 



(36) 
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rH2o = JAct-bz-»>r-bA(1.0 + P) 
ro2 = Daci G^ct + W2 ■ bz/2 - V2 



(37) 



-(1.0 + p)(a + bV2)-fe^^2 



+ mA,/2 + ONkBL(2(pR-1.0) 

S02 ^ (Tna + k^i) + (1.0 - 0|^)x^4^aj^4y^^ 

TnA ^ (^MAF-14 ^MAF-IS " ^MAFAe/^ ' ^MAF-6) 



(38) 
(39) 
(40) 
(41) 



In these relationships the subscript "Act'' means an effluent measurement or assumption (an 
"actual" value). The term J^^^^ in Eqs.(37) and (38) relating to the moles of effluent H2O could be 
input as a constant value or measured. All other values in Eqs.(36) through (41) are either 
evaluated explicitly based on input data, intemal models and/or have minor import but are carried 
in the formulations for consistency. The amount of inert matter (xttjo) associated with black liquor 
fuels as used in recovery boilers is typically small, less than 5% by w^eight and may be held 
constant which is the Preferred Embodiment, or as may be determined using the Air/Fuel ratio as 
taught in '994, or as otherwise obtained. 

[0052] As a group, these relationships are of critical importance for understanding The 

Input/Loss Method as applied to recovery boilers. If fuel chemistry is to be resolved, thus heating 
value, boiler efficiency and thus accurate fuel flow, etc., then stoichiometric relationships generally 
representing Eqs.(30) to (41) must be resolved. These equations are not unique in their grouping of 
terms, further reductions and/or complexities are possible. The grouping of terms adopted here 
principally follows from the right-side of Eq.(19BL). 

[0053] For the foUov^g discussion, assume initially that tube leakage, fuel inerts and fiiel 

inorganics (bp, a^Q, ^14, and a^^) have zero values. With this assumption, Eqs.(30) through 
(34) yield five equations with eight unknowns. For this situation, unknowns include through a^, 
and the terms "a" and "x'\ The term ''x'' is a convenience term and could be divided through 
changing the base of Eq.(19BL) to unity moles of fuel, thus eliminating use of xoj terms 
comprising two unknowns. However, if done, then the effluent's base becomes per mole of fuel, 
e.g., D^^./x, adding a different complexity involving the normalization of effluent measurements. 
Although the requirement £aMAF-j ^ ^ -^^ a convenience, it affords another, and viable, equation. 
By making a molar nitrogen balance, and assuming 100 moles of dry gaseous effluent at the 
boundary, the "a" quantity (moles of combustion oxygen) may be resolved independent of Eq.(30), 
thus reducing the unknowns; detailed below. Again, the entire combustion equation, Eq.(19BL), 
could be divided through by a^, or xa4, setting a carbon base. Effluent N2 could be resolved by 
difference assuming 100 moles of gaseous effluent (CO2, HjO, O2, SO2, the minor pollutants being 
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measured or assumed), or N2 could be measured directly. Or, further still, by assuming constant 
values for fuel nitrogen and sulfur, a| & (typically minor fuel constituents), with resolution of 
"a", and say: a3 = 1.0 - E^MAF-j J the system is reduced to three equations with four unknowns; 
these include Eqs.(31) through (33), with a^, and 'V. As another example, if aj4, 
and are assumed constant, then the combined Eqs.(31) and (32) (with cancellation of xa2) 
represents one equation with two unknowns, ''x" and however at least and a 14 represent 
major and variable constituents of black liquor fuel. And, of course, further reductions and 
manipulations of unknowns and equations is entirely possible. However, close examination of the 
physical problem of combustion stoichiometrics, in which fuel chemistry is to be determined from 
Choice Operating Parameters, indicates that the mathematical system invokes a stifiF matrix - and of 
course having fijrther complexity if the black liquor's inorganic components are fully involved. In 
summary, these manipulations are discussed to emphasize that, as taught by this invention, 
algebraic manipulations must address the physical reality of the thermal system, this means the 
fuel's intrinsic chemical relationships. Such intrinsic chemical relationships must recognize for 
example that black liquor is a highly wet fiiel consisting of both hydrocarbons and inorganics, and 
must also recognize the influences played by the fiiel on the smelt measurements (p^ and (p^ by the 
principal gaseous effluent measurements D^^^, J^^^ and Gjs^^. Invoking such intrinsic chemical 
relationships reduces The Input/Loss Method's sensitivity to instrumentation errors and enhances 
the accuracy of the applied minimization techniques which correct such measurements. The 
solution problem is to employ intrinsic chemical relationships within the resolution of the 
combustion equation. 

[0054] To address the solution problem, the Preferred Embodiment capitalizes on the 

intrinsic chemical relationships found in organic fuels between MAF hydrogen and MAF carbon, 
and, for black liquor fuel, relates stoichiometrically the organic to the inorganic compounds 
through the Sodium to Carbon ratio ((p^)- Further, the fiiel mole term, x, is used in the Preferred 
Embodiment for this single variable appears in all stoichiometric conversions to mass flows (the 
consistent determination and conservation of which is an object of this invention); and represents 
an iterative parameter for Eq.(19BL). Further, the Preferred Embodiment does not require that the 
minor fuel constituents be assumed constant. Further, as will become apparent, the Preferred 
Embodiment allows use of multidimensional minimization techniques which addresses 
instrumentation errors. 

[0055] Returning to the solution problem as posed by Eq.(19BL), if fiiel chemistry is to be 

resolved then fiindamental problems require solution - independent of algebraic manipulations. 
These problems include the following: 1) the black liquor's inorganic components; and 2) for any 
set of unknown fuel quantities, a:, there is always, at least, additional unknowns: ia the combustion 
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air term, "a"; in the fiiel mole term, 'V; and in the unknown tube leakage (b^) a subject of this 
invention. The first of these problems is resolved through the Sodium to Carbon ratio (pj^j, and 
obtaining correlations of molar MAF fractions of fiiel potassium and chlorine as a fiinction of 
sodium: 

^MAF-15 ^ ^15 Sl5^1VIAF-14 (42) 
^MAF-16 ^ + Bi5a]vjAF-14 (43) 

[0056] The second of these problems is solved, in part, by reducing Oj quantities to a molar 

MAF bases, eliminating the influence of the two components not chemically involved in the 
hydrocarbon nor inorganic portions of black liquor's make-up, water and inerts, and then solving 
for the remaining a^^.j terms. Dependent a^AF.j terms may be algebraically resolved by 
obtaining correlations of molar MAF fractions of fuel hydrogen, nitrogen and oxygen to molar 
MAF fixel carbon, Eqs.(48), (49) & (50), and molar MAF fuel potassium and chlorine to molar 
MAF fuel sodium, Eqs.(42) & (43); thus establishing intrinsic chemical relationships. Fuel carbon 
is then resolved explicitly by: converting the left-side of Eqs.(31), (32) & (33) to a MAF base; then 
subtracting Eq.(31) from (32) to eliminate the water term then converting inorganic 

compounds associated with T^a, using Eqs.(42) & (43) and the (p^ ratio, to MAF carbon terms; 
substituting Eqs.(49) & (50); resulting in two functional equations, the combined Eqs.(31) & (32) 
and Eq.(33), with two unknowns, Xj^ and a^AF^; solving produces: 

^MAF-4^ [" rco2(^C7 ^Cs) ^C4^C2] 

/ [W^cs + W ^ ^C4 (1-0 - ^C3)] (44) 

^MAF ^ ^C02 I [^MAF-4 " ^MAF-4^C3 ' ^cJ (45) 

However, before solving Eqs.(44) & (45), the "a" quantity (moles of combustion oxygen) used in 
and b^, and thus in ^^4, ^(-9 & Fj^o? resolved via Eq.(46) by recognizing: 1) the bases of 
Eq.(19BL) is that, at the boundary: 100 dry gaseous moles = D^^^ + Qp^^ + (h -H pacj)^^^) + e^^.^ + 
^Act ^n'^bl (10 - ^n)^maf^maf-65 is used to substitute for "h" in a nitrogen molar balance 
simplifying with V^^^ and 2) that: pa = 100(/JAct ' l-0)/[^Act (1-0 + (|)Act)l- 

a =(rN2 - XMApaMAF-lV^Act (46) 



In these equations the following terms are developed from algebraic manipulations and 
simplifications. Furthermore, these equations example the process of establishing intrinsic 
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chemical relationships as developed with the objective of establishing functionality between 
Choice Operating Parameters (e.g., effluent data) and fuel chemistry. 



^C2 ~ [^15 " ^16^2 

- ^sC^MAF-io + + Ai6/2 + v/xmaf)] (47A) 
5c3 = aN(PN[1.0 + B,5-Bi6/2 

-<Ps(1.0 + Bi5 + Bi6/2)] (47B) 

^C4 = - W2 + + 2ONkBL(1.0 - (Pr) 

+ (1 .0 - aN)xMAFaMAF-6 (47C) 
5c5 = aN[1.5(Ai5.Ai,/2) 

+ (0.5 - 2<|>R)(ps(aMAF.io + + Ai6/2 + v/x^af)] (47D) 
^C6 = aN[L5(pN(1.0 + Bi5-Bi6/2) 

+ (0.5 - 2(|>r)(Ps(|>n (1.0 + B,5 + Bi6/2)] (47E) 

^C7=-A3 + A5 (47F) 

^C8=-B3 + B5 (47G) 

+ ^MAF^MAF^ (<^N " ^OnCRr + 1 .0) (47H) 



[0057] Note that the v term of Eqs.(19BL) & (39), relating free carbon found in the smelt, 

may be assumed constant or may be expanded as a function of a^^^ or aMAF,i4 if its magnitude 
effects the accuracy of the computed aj^^^^. Given CLyy^_^, ^maf resolved via Eq.(45), 

based on Eq. (33). Given a^^j^^, sodium, aMAF-i4' then resolved from ((p^ <X'MAF-d' Note that a 
system carbon balance involves all major stoichiometric terms including x^^, ^Act^ ^Act 

and p. In Eqs.(42) & (43) the constants Aj and Bj may be developed as part of the Reference Fuel 
Characteristics associated with a specific black liquor taken from the actual As-Fired, or obtained 
from generic specification of a black liquor fiiel, or otherwise obtained. Such data may also be 
collected, or otherwise obtained, for fiiel nitrogen, oxygen and hydrogen. Note that if ultimate 
analyses are not available for the fiiel, for potential used in developing Aj and Bj, and especially A5 
and B5 for MAP hydrogen, then determination of elemental (ultimate) analyses may be had from 
so-called proximate data. Such conversions from proximate to ultimate analyses is well established 
art, four such conversion methods are described in A.B. Gill, Power Plant Performance. 
Butterworths: London, 1984, Chapter 2, pages 70-77. And although such proximate to ultimate 
conversions are not the Preferred Embodiment, it affords another example of establishing intrinsic 
chemical relationships required of Reference Fuel Characteristics. 
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aMAF-i=Ai+BiaMAF-4 (48) 

aMAF.3 = A3 + B3aMAF-4 (49) 
^MAF-S " A5 + BsttMAF^ (50) 

Fuel MAF sulftir is resolved explicitly by solving Eq.(34) and applying Eqs.(42) & (43) and the 
ratios <ps and (p^ : 

^MAF-6 ^ ^S02 I ^MAF (51) 
= (Ps [ttMAF-lO Ai5 + A16/2 + v/Xmaf] 

% <Pn aMAF-i [1-0 + Bi5 + + v/xmaf] 

+ aNkBL/XMAF (52) 

w^here the inorganic terms comprising T^q2^ whose compounds describing the inorganics 
produced from black liquor combustion, see Eq.(19BL), are determined as foUov^s: 

^n^ko ^maf ^maf-15 (53) 

%Tnc ^maf ^maf-16 (54) 

% T|^s 9r (^maf ^maf-6 ■ '^bl) (55) 

(1.0 - (Pr) (Xmaf ttMAF^ ■ ^bl) (56) 

%(TnO'^1^Bl) ^^MAF^na ^n'^BL ■ %^MAF^MAF-15 (57) 

After Eqs.(44), (49) & (50) are resolved, fiiel MAF water is resolved explicitly by employing 
Eq.(3 1), or by adding Eqs.(3 1) & (32) for numerical over-check, through substitution of Xj^^: 

^MAF-2 ^ (PmO ' ^MAF ^MAF-5) / ^MAF (58A) 
^ (" ^MAF-5 ^9 ^MAF-S^* H2o) ^ (^C9 " ^md^) (58B) 

[0058] Fuel MAF nitrogen concentration is small, and typically may be fixed as a constant. 

Fuel oxygen may be described by Eq.(49). The MAF sulfiir concentration may be computed 
assuming: aMAF-6 ^ ^6^uafa^ ^ opposed to Eq.(52), if the effluent SO2 measurement is 
questionable (or if limestone is injected and its conversion rate is questionable). As another 
altemative, either the oxygen or sulfur could be used to assure that Eaj^i^.j = 1.0, where: j = 
1,3,4,5,6,14,15,16. The sole criteria in deciding the exact methodology is the reliability and 
availability of effluent data and its relative impact on MAF fiiel terms. The Preferred Embodiment 
is to solve OLyy^_^ using Eq.(49), a,^^^^ using Eq.(52), and a^^.j by balance. For example. 
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although nitrogen is a major combustion effluent, using Eq.(30) to solve for fiiel nitrogen, typically 
a very minor component of the fueFs makeup, would invite even slight errors made in determining 
effluent N2 (made either by direct measurement or determined by difference suggested in Eq.(36)) 
will greatly amplify the uncertainty in aMAF-i- Indeed, the term is used to resolve the 
combustion O2 term via Eq. (46), wherein the term Xj^y^aj^i^.j is resolved by balance and iteration 
on Xj^. Besides carbon, an important term is the fiiel MAP hydrogen concentration which 
represents an intrinsic chemical relationship with carbon. Experience has shown that a valid 
relationship may be achieved since the hydrogenxarbon chemical bond is predominant with 
carbonxarbon. With an established MAP hydrogen to carbon relationship, the variability of a 
specific black liquor fuel then lies with its water, inorganics and inerts. 

[0059] To this point in the process, resolution of fuel chemistry has been possible by 

assuming a Moisture-Ash-Free (MAP) basis. In carrying this out, it is not critical to the invention 
that carbon be computed as the independent quantity; nor is it critical to use a MAP basis as the 
inert fi'action is typically small in black liquors, thus a dry basis may serve equally well. Hydrogen, 
or any other MAP fiiel element could be considered as the independent variable. Por example, if 
solving for hydrogen first, then a^^_^ =/ (aMAp-s)- However, the Preferred Embodiment is to 
place independence with carbon, where of course the greatest accuracy and sensitivity may be 
found in the effluent measurements. CO2 does not exist in the combustion air to any appreciable 
concentration; it does not leak into the system; it is generated only fi-om combustion and as a major 
effluent has obvious sensitivity to resolution of a|yi^^ of Eq.(44). 

[0060] Thus all fuel constituents are therefore determined on a MAP bases. Prom these 

values, the wet base molar fuel fractions are then determined, as are the wet base moles of fuel (x) 
and the wet base (As-Pired) weight fractions (WF^ of all fuel components j: 



a j = ^MAF-j / (1 -0 + aMAF-2 + aMAF-io) (59) 

X = Xmaf (10 + aMAF-2 + ocmaf-io) (60) 

^^j ^ ^MAF^MAF-j (61) 

FFFj=ajNj/(EajNj) (62) 

^DRY-j = WF ^/ (1.0 (63) 



[0061] As taught in '879 the sensitivity of computed fuel chemistries to an effluent H2O 

measurement is substantial When employing any of the Input/Loss methods which compute fixel 
chemistry, fiiel chemistry will effect the computed fiiel heating value, boiler efficiency and thus 
system eflBciency. With such Input/Loss methods, computed fiiel chemistries are generally 
sensitive to combustion gaseous eflQuents, including CO2, H2O and O2. The effluent H2O 
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measurement (or its assumption), as well as effluent CO2 and O2, may be corrected through 
methods taught by '879; however, such measurements will be impacted by tube leaks, and tube 
leaks of even moderate flows. Both '879 and this invention discuss the optimization of certain 
parameters, such parameters include an assumed flow of water in-leakage into and mixing with the 
combustion gases; thus the modeling of tube leaks using a combination of combustion 
stoichiometrics, Eq,(19BL), etc., and optimization methods. 

[0062] Having computed black liquor's fuel chemistry, the heating value of black liquor is 

next computed. Following the teachings of '994, black liquor heating value is determined base on 
a differential analysis. Note that the term Nj^^^^ is the molecular weight of the MAF-base fuel 
(without water and inerts). 

AHHVmaf^i^ = HHVMAF-Ref - (- 178387.18aMAF-3 + 183591.92aMAF^ 



+ 78143.68aMAF-5 " AH^^al/BLW / NMAF-Ref (64) 

HHVMAF-uncorr = (" 178387.1 8aMAF.3 + 183 591. 92aMAF4 

+ 78 143.68ai^„5 - AH^^al/BOActual / NMAF-Actual (65) 
^Hc-Cal/BL ^ [NNa2S04^J^-Caima2S04 

+ Npj2sAH^^ajyj^S -2NNaHS^H/^<:aJ/NaHsl (66) 

HHVmaf = + AHHV MAF-delta (67) 

HHVdry HHVmaf(10 - WF^KYAo) (68) 

HHV^ =HHVDRY(10-f^2) (69) 



Eq.(66) illustrates, for example, that black liquor may contain only Na2C03 and K2CO3 (which are 
assumed to appear as products),and NaHS which reacts to produce Na2S04 and H2S, thus requiring 
a AHeat of Combustion correction, AH^^^j/bl- In Eq.(66) MAF sulfiir is first taken as the 
reference value, oc^|y^_5_Rgf5, wherein AH^^^l is then used in Eq.(64); then it is taken as the 
actual, a^AF^ via Eq.(52), wherein ^H^^qUbl is then used in Eq.(65). 

Boiler Efficiency for Recovery Boilers 

[0063] Boiler eflSciency computations are a vital portion of any of the Input/Loss methods, 

and especially so when resolving tube failure leakage flow. Methods presented in '429 describe 
how to compute a high accuracy boiler efficiency applicable for conventional steam generators 
burning common hydrocarbon fuels. However, as applied to recovery boilers, given black liquor 
fuel contains reactive inorganic compounds subjected to reduction, the methods of '429 require 
modification. Methods presented herein describe how to compute a high accuracy boiler efficiency 
applicable for black liquor fuels. 
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[0064] Notably The Input/Loss Method emphasizes computational consistency, consistency 

principally achieved through boiler eflBciency and its computation of heating value as dependent on 
the calorimetric temperature, Tq^. The calorimetric temperature must be consistently employed in 
all energy terms and Heats of Formation. For recovery boilers is it common industrial practice to 
correct the measured heating value for Heats of Formation associated with reactions particular to 
recovery boilers, and commonly the chemical reduction of Na2S04. Such industrial corrections are 
taught in the following documents: "Performance Test Procedure Sodium Based Recovery Units", 
CA Report No. 84041601, March 1996, TAPPI Press, Atlanta Georgia (see paragraphs 5.2.2.5, 
7.2.2.8 and 7.2.2.9, and Appendk B); and Steam, Its Generation and Use , 40th Edition, Edited by 
S.C. Stultz and J.B. Kitto, published 1992 by The Babcock & Wilcox Company, Barberton, Ohio 
(see Chapter 26, page 26-5). Such common industrial corrections address the difference between 
ideal combustion products associated with a bomb calorimeter versus actual products associated 
with further reducing black liquor compounds. However, as taught herein such industrial 
corrections are thermodynamically inconsistent, and simply wrong. Industrial corrections reduce 
the higher heating value with a computed AH^ term: HHVP - AH^; this term effecting the 
denominator in conmion boiler eflBciency formulations. Given the teachings of '429 the AH^ term 
is an intrinsic portion of the numerator of boiler efficiency; an affect of describing ideal versus 
actual products of combustion bound in the HPR^^^ and HRX^^^^ terms. Further, the principles 
taught about the use and computation of L Factors, an important System Effect Parameter concept, 
are fundamentally based on the as-measured heating value combusted under theoretical conditions 
with dried foel (not a HHVP - AHr concept); see the section on "System Effect Parameters" and its 
Eqs.(272) & (273) below. '429 teaches, for a HHV-based eflBciency, that: 

%-HHV = BBTC 

mAF(HHVP + HBC) 

= HHVP + HBC - ^Losses/m >^ ^j22) 
HHVP + HBC 

^ HHVP + HBC - (HSL + HNSU ^^23) 
HHVP + HBC 

where the numerator of Eq.(123), if considering only HSL losses for clarity (i.e., ignoring 
Non-Chemistry & Sensible Heat Losses, 11^=1.0), may be expanded as: 
= HHVP + HBC-HSL 

= - HPRA^jt.HHY + HRXy5^^^.HHV 
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= - HPRaci.hhv HRXcaj.HHV + 

= -HPRAct-HHV + HPRMeal-HHV + HHVP + HBC 

where this last expression illustrates that actual product compounds (e.g., as eflfected by the 
process' Reduction Efficiency) are net of ideal combustion products (found in a bomb calorimeter 
producing the measured HHVy^). Heats of Formation, and sensible heats associated with the 
Firing Correction term HBC, are all thermodynamically consistent When reducing Na2S04 in a 
recovery boiler via the reactions: 

Na2S04 -> Na2S + 2O2 (124) 
CO2 + Na2S04 -> SO2 + Na2C03 + O2/2 (125) 

the Heats of Formation of the non-ideal product compounds Na2S, SO2 and Na2C03 are addressed 
through HPR^ct and its Eq.(133), whereas the ideal - and the measurement base for all heating 
values - are addressed through HPRy^^. Note that Na2C03 as produced from reducing Na2S04 is 
govemed by the effluent measurement of SO2 O^bl)^ although may be inaccurately measured, 
the measurement may be corrected or otherwise obtained. The quantity BBTC/(HHVP + HBC) of 
Eq.(121) must be constant for a given in-situ system burning a defined fuel; thus any changes in its 
combustion process (e.g., changes in Reduction Efficiency, <pj^) which will of course affect fuel 
flow, must also affect t|b (through HPR^^^ in proportion, mamtaining BBTC/(HHVP + HBC) 
constant. A (HHVP - AHr) procedure is inconsistent since the term: BBTC/(HHVP - AHr + 
HBC) would not remain constant since AHr =/ ((p^); the term would produce a variable fiiel flow. 
[0065] Non-Chemistry and Sensible Heat Loss term (HNSL) is computed as taught in '429. 

With its determination, only the three major terms HPRA^.t, HRX^^.^ & HBC remain to be defined to 
complete boiler efficiency. These are taught in the following paragraphs. To fiilly understand the 
formulations comprising HPR^^^ HRX^^^^ and HBC, take note of the subscripts associated with the 
individual terms. As example, when considering water product created from combustion, nQjmb-mo 
of Eq.(131), its Heat of Formation (saturated liquid phase) at T^ai must be corrected for boundary 
(Stack) conditions, thus, hst^ck " ^^.^^1- Th^ Enthalpies of Reactants of Eqs.(135) & (136) are 
determined from ideal products at T^^, the Firing Correction then applied. 
[0066] Differences in formulations required for higher or lower heating values should also 

be carefully reviewed, as recovery boilers are evaluated using either. Higher heating values require 
use of the saturated liquid enthalpy evaluated at T^^^; lower heating values require the use of the 
saturated vapor at Tq^, The quantities which are not so corrected are the last three terms in 
Eqs.(13 1) & (132): water bom by air; known in-leakages (b^); and tube leakage (b^). These terms 
undergo no transformations since they have non-fiiel origins. Heating values and energies used in 
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Eqs.(131) through (137) are always associated with the system boundary, specifically defined by: 
the entry point of the As-Fired fiiel (or the "supplied" fiiel in the case of fuel rejects); the ambient 
air conditions; the Stack location including gaseous effluent temperature and the Continuous 
Emission Monitoring System (typically measuring CO2). 

[0067] Boiler efficiency is defmed as either HHV- or LHV-based, being composed of a 

combustion efficiency and a boiler absorption efficiency: 

TlA (126) 
%.LHV = 'nC-LHvTlA (127) 

The boiler absorption is determined from the Non-Chemistry & Sensible Heat Losses term (HBC) 
as fully discussed in '429. It is an important teaching herein that fixel flow compute identically 
from either efficiency base; thus, assuming no tube leakage: 

°*AF = TiB.HHN^VP + HBC) = riB-LHv(LHVP + HBC) (128) 

Such computations of fiiel flow using either efficiency, at a defined Tq^, is an important numerical 
overcheck of this invention. 

[0068] Combustion efficiency is deteraiined by the following, as either a HHV-based or a 

LHV-based efficiency: 

The development of the combustion efficiency term, as computed based on HPR^^^t & HRXy^^t and 
involving systematic use of a combustion equation, such as Eq.(19BL), is believed an improved 
approach versus the primary use of individual "stack loss" terms. Mis-application of terms is 
greatly reduced, while numerical accuracy is increased Most importantly, valid system mass and 
energy balances are assured. 

[0069] The Enthalpy of Products (HPRy^^t) term is as follows. For higher heating value 

calculations: 



[ ^mb-H20 (^I^-Cal/liq *^Stack " *%Cal) 
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%uel-H20 (^Stack hf^al) 
^Air-TOO C^Stack ■ '^-Cal) 
(hstack ■ hsteamz) 

+ W (hst3^i, - hstea^j) ]h2o Nh2o /(xNaf) 



For lower heating value calculations: 



^ [ ^omb-H20 (^H/-Cal/vap ^stack 
^ %uel-H20 (^Stack " hg^al) 
^Air-H20 (^Stack ■ hg-Cal) 
(^Stack " l^Steamz) 



^-Cal) 



+ W (hstack - hsteamT) 1h20 Nh20 'i^P^, 



w4iere: HPR -^ = Enthalpy of non-water product i at the boundary ; 

refer to the right side of Eq.(19BL); for 
example, if product i is Na2S then: 

'^Na2S ^ Tns ^ 9r (^MAF ^MAF-6 ' ^^Bl)- 



HPR , . [AHj^3,i + Jdhi]niN/(xN^) 



^ stadc 



'^omb-HZQ 
%iiel-H20 

^Air.H20 



^Stack- 



H20 



**SteamZ 
l^SteamT 



^Cal 

= Molar water found at the boundary formed directly 
from combustion of a hydrocarbon (xa5). 

= Molar water found at the boundary bom by As-Fired 
fuel (as total inherent and surface moisture). 

^JAct-[bA(l-0 + P) + bz + br + xa5]=xa2 

= Molar water found at the boundary bom by 
combustion air and air in-leakage. 

=/(Pstack-H20^Tstack). where Pstack-H20 water's partial 
pressure per total wet molar and is given by: 

Pstack-ffiO = PAmb(JAct ^ PbJ/ClOO + J^^ + PbJ. 

= Weighted specific enthalpy of b^ in-leakage. 

= Specific enthalpy of tube leakage (br). 

= Heat of Formation of compound i evaluated 
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at the calorimetric temperature, Tq^. 

.T-Cal /.T-Cal 



dhcompound/i ^ L J ^U^^s (134) 
77 77 

[0070] The Enthalpy of Reactants (HRX^ct) term is as follows. Note that although SOj is a 
common gaseous product of ideal combustion of conventional fossil fuels, it is not assumed to be 
produced from the ideal combustion of black liquor. For higher heating value calculations: 

HRXAct-HHV ^ HHVP + HBC 

+ HPRcOl-ldeal ^^^H20-IdeaI-HHV 
+ iZP/?l^a2C03-Ideal ^^^Na2S04-Ideal 

+ ffi^^NaCl-Ideal ^^^K2C03-Ideal (135) 

For lower heating value calculations: 
HRXAct-LHV^LHVP + HBC 

+ HPRQQ2-UeQ\ ^^^H20.Ideal-LHV 
+ /rP/?Na2C03-Ideal ^ ^^^Na2S04-Ideal 

+ ^Pi^NaCWdeal ^^^K2C03-Ideal (136) 

where: ffPRco2.\dsai ^ Energy of CO2 ideal product from complete 

combustion at the calibration temperature. 

^^^H20-ideai-HHV Energy of H2O ideal product from complete combustion, 

ending with condensed water, at the calibration temperature. 

^(^-Cal/liqN)H20a5/NAF 

//P/?jj20.ideai-LHV " Energy of H2O ideal product fix>m complete combustion, 

ending with water vapor, at the calibration temperature. 

^(^-Cal/vapN)H20a5/NAF 

^^^Na2C03-ideai ^ Energy of Na2C03 ideal product from complete 

combustion at the calibration temperature. 

= ^f^-Cal/Na2C03 ^ T^q N|^a2C03 ^ (^Af) 

^^^Na2S04-ideai ^ Energy of Na2S04 ideal product from complete 

combustion at the calibration temperature. 
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- ^H/.Cal/Na2S04 ^nTna NNa2S04 ^ (^Af) 

^^^Naci-ideai ^ Energy of NaCl ideal product from complete 

combustion at the calibration temperature. 

= AH^.C3,/NaCl %Tnc N^aCl / (?^Af) 

^^^K2C03-ideai Energy of K2CO3 ideal product from complete 

combustion at the calibration temperature. 

= AI^^3j/K2C03 TkO Nk2C03 ^ (^Af) 

It should be noted that in Paragraphs 0069 and 0070 when using the fuel's calorimetric temperature, 
established when determining the ftiel's heating value, as the thermotfynamic reference energy 
level for tiie Enthalpy of Products and for the Enthalpy of Reactants, that these include corrected 
Heats of Formation, AH^^^j/i for example as taught by Eq.(134), said Al^^^i/i term being used in 
Eqs.(131), (132), (133), (135) and (136). Further, that AH^^ai/i is corrected relative to a standard 
Heat of Formation for substance i, taken at 77 F (25 C) and denoted by AH^.^^/j . 
[0071] The Firing Correction (HBC) term is given by the following. 

HBC^ Cp(Taf - Tcai)Fuel + (QsAH + WpD)/mAF 

+ [(hAmb - hcal)Aira(1.0 + P)(1.0 + (t)Act)NAir 
+ (hg-Amb - hg<:al)H20 ^aCI -0 + P)Nh20 

0^SteamZ"hf-Cal)H2o'^zNH2O 
+ (hsteamT " hf<:al)H2ol>rNH20 1/(xNaf) (137) 

where: hg.Amb-H20 ^ Saturated water vapor enthalpy at ambient diy bulb, T^mb- 
Cp(Taf - Tcai)Fuei Sensible heat in the As-Fired fiiel relative to Tq^. 

(^Amb " ^ai)Air ^ AEnthalpy of combustion diy air relative to Tq^. 
0^-Amb ■ '^-Cai)H20 ^ AEnthalpy of moisture in combustion air 

relative to saturated water vapor at Tq^. 
O^steamz ' ^^f-cai)H20 ^ AEnthalpy of bz known water in-leakages 

(at an average h^teamz) relative to 

saturated liquid water at Tq^, 
(hsteamT ' ^f^hio = AEnthalpy of tube leakage (at hsteamx) 

relative to saturated liquid water at Tq^. 



[0072] The above equations are dependent on common system parameters. Common 

system parameters are defined following their respective equations, Eqs.(131) through (137). The 
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BBTC term, comprising common system parameters, is determined from commonly measured or 
determined working fluid mass flow rates, pressures and temperatures (or qualities). Further, 
supporting terms such as thermocfynamic properties, radiation & convection loss curves, guidelines 
for estimating Non-Chemistry & Sensible Heat Losses, sodium compound effects, miscellaneous 
terms, etc. are discussed in the following codes, standards and patents: the American Society of 
Mechanical Engineers' (ASME) Performance Test Codes (PTC) 4.1 and 4; the "Performance Test 
Procedure Sodium Based Recovery Units", CA Report No. 84041601, March 1996, TAPPI Press, 
Atlanta Georgia; the German standard "Acceptance Testing of Steam Generators, DIN 1942, DIN 
DEUTSCHES Institut Fur Noraiung E.V., February 1994; the European standard (draft) prEN 
12952-15:1999 (also: CEN/TC 269/WG 3 N 337), "Water-Tube Boilers and Auxiliary Installations 
- Part 15: Acceptance Tests", November 1999, European Committee for Standardization, Central 
Secretariat, rue de Stassart, 36, Brussels; the British Standard "Code for Acceptance Tests on 
Stationary Steam Generators of the Power Station Type", BS 2885:1974, ISBN: 0 580 08136 2; 
and throughout US Patents '420, '470, '994 and '429. If conflicts arise between any method or 
procedure of these codes, standards or patents, this disclosure's methods, procedures and invention 
shall prev£uL 

[0073] It is interesting that the German standard DIN 1942 states that its reference 
temperature (t^) for boiler efficiency is 25 C (77 F). However, in its paragraph 6.2, DIN 1942 
allows that "other temperatures may be agreed upon" specifying in its Eq.(lb) how to correct 
heating value for a (t^ - 25) effect. DIN 1942's Eq.(lb) not only corrects heatmg value using 
sensible heats (CpAT), incorporating all terms described by the Firing Correction (HBC) taught 
herein, but also corrects heating value for flue gas sensible heat and thus the denominator of rig. 
This flue gas correction is incorrect. As taught herein a flue gas sensible heat term appears in 
Eq.(133), integrating the (dh) energy of combustion products from Tc^i to T^^^i^ effecting the term 
(-HPR^jjt.LHV HR^Act-LHv)^ 3nd thus only the numerator of t|b. Further, DIN 1942 does not 
teach to correct the Al^^^i/i term of Eq.(133) for a reference temperature different from T^ai or t^,; 
the term is ignored in DIN 1942 but is taught herein through Eq.(134). DIN 1942 does not correct 
Heats of Formation to achieve consistency with calorimetric temperature as taught herein; indeed, 
nowhere in DIN 1942 does any Heat of Formation appear. 

[0074] The draft European standard mimics DIN 1942; it also discusses "Input/Loss 
Calculations" in its Section 9.4.5.2, but these calculations are not related to "Input/Loss methods" 
as used herein as defined in paragraph 0035. The British Standard requires that 'the standard 
reference temperature for the calorific value of solid, liquid and gaseous fiiels shall be 25 C" 
(Clause 2), and offers no means to alter. The TAPPI performance test standard references energy 
flows to a constant 77 F, but offers no means to alter. TAPPFs Appendbc B develops a "Heat of 
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Reaction Correction", correcting heating values as (HHVP - AH^^), however no T^ai correction or 



PTC 4.1 references energy flows to an arbitrary ''reference air temperature", T^^. PTC 4 
references energy flows to a constant 77 F, but offers no means to alter. Neither DIN 1942, the 
draft European standard, the British standard, the TAPPI standard, PTC 4.1 nor PTC 4 mention 
how a reference temperature should be established, other than setting it constant. There is no 
mention of the fiiel's calorimetric temperature in any of these standards, nor the correction of 
reference Heats of Formation to the calorimetric temperature, nor using the fuel's calorimetric 
temperature as the reference temperature. Note that Heats of Formation are normally cited in the 
literature as being referenced to 25 C and 1.0 bar pressure. The Preferred Embodiment of '429 and 
for this disclosure requires the reference temperature to be taken at the calorimetric temperature 
established when determining heating value, and that all energy terms, including those associated 
with Heats of Formation, HPR^j^^t, HRX^^^ and HBC, be consistently referenced to this same 
calorimetric temperature. 

[0075] Several PTCs and "coal" textbooks employ simplifying assumptions regarding the 
conversion of heating values. For example, a constant is sometimes used to convert from a 
constant volume process HHV^ (i.e., bomb calorimeter), to a constant pressure process HHVP. 
The following is preferred for completeness, for solid and liquid fuels: 



where, in US Customary Units: T^aj^bs absolute calorimetric temperature (deg-R); Rq3s = 
1545.325 ft-lbfi^mole-R; and J^ech = 778.169 ft-lbfi®tu. For gaseous fiiels, the only needed 
correction is the compressibility factor (2) assuming ideally computed heating values at the 
specified Tq^: 



Z and HHVjdeai may be evaluated using American Gas Association procedures. To convert from a 
higher heating value (also termed gross or upper) to a lower heating value (also termed net) use of 
Eq.(142) is exact, where ^fg^jQaumo is evaluated at Tq^. The o>o^gen in the effluent water is 
assumed to derive from combustion air and not from fuel oxygen: thus is not included as there is 
no molar change, fiiel oxygen is taken as the diatomic in Eq.(19BL). 




HHVP =HHVaf + AHv/p 

^V/P - Rcas TcaUbs (^5^^ " «l) / (5reJ ^Mech ^af) 



(138) 

(139) 




(140) 
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LHV^=HHV^ - AHi^ (141) 

AHl^ = ^fg^al/H20 ("2 + a5)NH20/(^REJ Naf) (142) 

Within Eq.(142) the ^^gj term accounting for rejected fiiel is defined as: = (1.0 - J^asH-af) 
/(LO - ^Ash-Sup)- ^REJ ^so corrects Eq.(144). These same procedures are applicable for a fuel 
cleaning process where the fuel's mineral matter (inert) is removed. 

LHVP =LHVaf + AHv/p - AH^^.^hv (143) 

AHeo„.LHV = AHi^(Crej - 1.0) / ^REJ (144) 

Tube Failure Detection Methods 

[0076] This invention teaches to add, apart fi-om b^ a term descriptive of tube leakage to 

the typical combustion equation; its symbol herein is b^ whose units are moles of liquid water or 
steam m-leakage per 100 moles of dry gaseous effluent, A procedure is taught by which the flow 
of in-leakage associated with a tube failure and its location may be determined This procedure 
relies on a hydrogen stoichiometric balance incorporating bp, applying limit tests on fuel 
constituents and tests on the rate of change of certain parameters; which tests provide an indication 
of possible tube failure and its stoichiometric causality (termed a "Tube Failure Mechanism", 
indicated by a numerical identification). Once a Tube Failure Mechanism has been identified, a 
tube failure flow rate is determined, but determined using a separative analysis technique. In 
summary, the Preferred Embodiment relies on modification of a traditional combustion equation, 
employing the hj term, it relies on limits and rates testing, and it relies on successive optimizations 
of Choice Operating Parameters which minimize errors in System Effect Parameters. All of these 
topics are discussed in detail in this disclosure. Forming a hydrogen stoichiometric balance using 
Eq.(19BL), and solving for hj, results in Eq.(200) and thus allows the determination of tube 
leakage in moles. Note that the effluent moisture (at the Stack), is defined as: Ja^ - j + PbA- Also 
note that Eq.(19BL) represents a mathematical model of the combustion process using a molar 
base. The manipulation of such a combustion equation when employing traditional fossil fuels 
(non-black liquor), is taught in '994; whereas methods required for black liquor fuels, considered 
novel, are as taught herein. 

br = Jact - x(a2 + a5) - b^ - bA(1.0 + P) (200) 



Eq.(200) illustrates that for hj to be positive, i.e., a tube leak being detected, that unique balance 
must be developed between the assumed (or measured) effluent water (J^^ and the predominating 
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negative terms: combustion water (xa2 + xa5), b^, and moisture in the combustion air and in the air 
leakage b^(1.0 + P). Eq.(200) demonstrates that use of an effluent H2O instrument, measuring J^^t* 
may not detect tube failures. For example, any unusual increase in J^^ could be caused by off- 
setting effects from high fuel water, high moisture in the combustion air, high £ur pre-heater 
leakage (a high p) and/or periodic soot blowing flow and/or use of atomizing steam (b^). Further, a 
tube leak could exist when the J^^^ term is decreasing as caused, for example, by a large decrease in 
fuel water (when, at the same time, by is increasing). To resolve such difficulties, this invention 
teaches the use of Eq.(200) in conjunction with one of the Input/Loss methods in which fiiel 
chemistry, the a j terms, are determined. 

[0077] When Input/Loss methods compute fiiel chemistiy associated with black liquor, 
such chemistry will include at least the determination of fuel elemental carbon (a4), fiiel elemental 
hydrogen (aj), fuel water (a2% and fiiel elemental sodium (ai^). Typical Input/Loss methods will 
determine such quantities, in part, based on Operating Parameters including principal effluent 
concentrations (CO2, 02 and H2O), combustion air psychrometrics (leading to b/^), and any known 
water and steam flows such as soot blowing and atomizing of fiiel (h^. Further, The Input/Loss 
Method teaches when applied to recovery boilers, as a portion of Reference Fuel Characteristics, is 
to obtain intrinsic chemical relationships between fiiel MAF hydrogen and carbon, and fiiel sodium 
and carbon, for example for hydrogen: Wmaf-s " A5 + Bsa^AF^; a^^^^ " A5 + B5aj^^ + 
C5(aMAF-4)^; or in general: a^AF-s =/(aMAF^); example for sodium: aMAF-i4 = <Pn «maf-4; 

or m general: aMAF-i4 ^/ (^maf-4)- Such fimctional relationships have been found to be most 
useful for fossil fuels including black liquor fiiels. It has been found that the constants A5, B5, C5 
and <Pn are typically remarkably constant for a specific fiiel. Such correlations then establish inter- 
dependency between fiiel carbon, fiiel sodium, and all principal effluents, and thus, through 
Eq.(200) and resolution of Eq.(19BL), between the effluent concentrations CO2, O2 and H2O, and 
the important hj term. In this context resolution of inorganic stoichiometrics (i.e., sodium) is 
vitally important to resolution of the organics, and thus to resolution of the tube failure (hj) term. 
Given such inter-dependencies, it is most likely that when assuming hj = 0.0, when in fact a tube is 
leaking, one or more fuel molar quantities (ttj, 02, ... a^, ai5, ai6) will compute outside 
reasonability limits or even as negative values. E?q>erience applying this invention has taught that 
fuel water will commonly compute as negative, even with a moderate leak when initially assuming 
br=0.0inEq.(200). 

[0078] In like manner, and especially for small leaks (when assuming bj = 0.0), the fiiel 
carbon, hydrogen and sodium terms could exceed reasonability limits; where for example, 
assuming that the constant B5 is negative and C5 is zero (which is typical for black liquor fiiels): 
«MAF-5 < aMAF-s/min a^^^/or tt^AF^ > aMAF^/max- Such bchavior when usmg Eq.(200), when first 
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assuming bj = 0.0 and then evaluating for reasonability limits, leads directly to an indication of 
tube leakage. This process is termed a possible Tube Failure Mechanism, that is an indication of 
possible tube leakage has been found by applying stoichiometric considerations (min/max checks); 
however further processing is called for to determine its validity and, if a valid leak, then to 
determine its mass flow rate and the location of the leak. 

[0079] To fiilly expand the concept of Tube Failure Mechanisms, TABLE lA and TABLE 

IB teach 35 such mechanisms, identified by a number, of how a tube failure may be detected 
through stoichiometrics, knowledge of such detection mechanism being important to the system 
operator. TABLE lA presents static trip mechanisms, while TABLE IB presents dynamic trip 
mechanisms related to rates of change. Experience in demonstrating this invention at two large 
power plants has indicated that making assumptions as to "apparently" impossible Tube Failure 
Mechanisms is not advised. Thus both minimum and maximum trip mechanisms are demonstrated 
in TABLE lA, as all are tested when monitoring a thermal system on-line when invoking the 
methods of this invention (i.e., invoking The Input/Loss Method's "Tube Failure Model"). For 
example, a cursory evaluation would suggest that a low fuel water concentration (a^AF ^ 
WF^q) could not indicate a tube failure given the mechanics of Eq.(200). However, if the thermal 
system experiences a small but steadily increasing tube leakage The Input/Loss Method of '879 
could steadily correct effluent water concentration downwards, causing failure mechanism #41 or 
#51; or water correction factors might exceed a lower bound causing failure mechanism ID #71. 
But also, unplanned scenarios of how The Input/Loss Method is correcting effluent water and other 
Choice Operating Parameters could create unexpected Tube Failure Mechanisms through complex 
stoichiometric relationships. Such considerations thus call for a blanket examination of all trip 
mechanisms, which is the Preferred Embodiment. TABLE lA encompasses the most likely of 
stoichiometric mechanisms based on the experiences gained demonstrating this invention, while 
TABLE IB encompasses likely rates of change of relevant parameters. Note that several of the 
rate mechanisms indicated in TABLE IB were found by testing the methods of this invention at a 
power plant; specifically Tube Failure Mechanism #64, #121 and #131 were found unusually 
sensitive. However, in extending the teachings of this invention as taught through TABLE 1 A and 
TABLE IB, there are other Tube Failure Mechanisms which may become apparent without 
departing fi-om the scope and spirit of the present invention. 
/// 
/// 
/// 
/// 
/// 
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TABLE 1 A: Static Tube Failure Mechanisms 



Mech. 
ID 


Trip Mechanism 


Comments 


11 


^Act ^AcJ/min 


Effluent HjO concentration at Stack; min. Jy^^t likely. 


12 


•^Act ^ •'^Act/max 


Effluent HjO concentration at Stack. 


21 


^MAF-4 ^ °'MAF-4/min 


MAF molar fraction of carbon; min. aMAF-4 likely. 


22 


"mAF-4 ^ <*MAF-4/iiiax 


MAF molar fraction of fuel carbon. 


23 


Negative square root 


Resolution of MAF carbon (aj^^^^) may require solving a 
second order equation, thus the possibility of tnppmg on a 

IlCgaLIVC dL[UarC lUUL \^Ur dliy IlidLllCIIlaLlWal UlUClCnillllcUC^. 


3 1 


/V <1 /V 

^MAF-5 ^ ^MAF-S/min 


lYl/TLr lilUlar iXaCLlUIl Ui iUCl iiyuiugcii. 


32 


^MAF-5 ^ ^MAF-S/max 


MAF molar fraction of hydrogen; max. aMAF-5 likely. 


41 


^MAF-2 ^ ^MAF-2/inin 


MAF molar fraction of fuel water. 


42 


^MAF-2 ^MAF-2/max 


MAF molar fraction of fixel water. 


51 


^^H20 ^ ^^H20/min 


As-Fired weight fraction of fuel water. 


52 


^H20 > f^H20/max 


As-Fired weight fraction of fiiel water. 


61 


^MAF-1 ^ ^MAF-l/min 


MAF molar fraction of nitrogen in fixel; nitrogen is 
computed by balance: aMAF-i = 1-0 - LO-yu^.\^ k^i,2,io- 


62 


^MAF-1 ^MAF-l/max 


MAF molar fraction of mtrogen m fixel; mtrogen is 
compuiea oy oaiance. clmafa ^-^ L^MAF-k,k^i,2,io- 


71 
/ 1 


^ ^2S/min 


i^orrecuon lacior lor eniueni rX2\J* 


79 


^ *-2S/max 


\^%jjTCK^ii\jMi idiLLur iur cuiucnL n2v/. 


81 


^ ^IS/min 


Correction factor for effluent CO2. 


82 


^ ^IS/max 


Correction factor for effluent CO2. 


91 


^MAF-14 ^ ^MAF-14/inin 


MAF molar fraction of sodium in fuel. 


92 


^MAF-14 ^ ^MAF-l4/max 


MAF molar fiaction of sodium in fixel. 


99 


my > 0.0 


Tube leakage found in previous monitoring cycle. 



/// 
/// 
/// 
/// 
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TABLE IB: Dynamic Tube Failure Mechanisms 



Mech. 


X 1 iU IVlCl^XiCUilMll 


V^UIIIIIICI lib 


14 




T^fiti^ nf* rVmnof^ in utii^rttTAptAH A'FFIii^nt l-T C\ 
jvduv duuigv ill uiivurrcL/LCU ciiiuciii xx^w. 


IS 


MC T VHt > T imit 


ixoic ui cndiige in LorreLicQ ciiiuciu xi2v^* 


94 


Hrn VHt > T imit 


IVaLC Ui Llldllgc 1X1 UIILUITCCLCU ClllUCIll \^W2» 




d(r T% VHt > T imit 
V IS Act/ IjUIIII 


jvaic oi Ciiongc in coiTccicci ciuucni v^v^2* 




A(C\ VHt > T imit 


ivdic OI cndngc in uncoircciea cmucni 


f\f\ 

uu 


A(C C\ VHt > T imit 


iVdic oi cndnge in cuircL'Leu cmucni yJ2* 




d(C \ld\ > T imit 


ivalC Ui L/IUulgC ill UlC C/UiTcL'UOn lO ClliUwni ^^2' 


lA 


H/T* ^/Ht > T imit 


ivaie OI cnougc in me correciion lo emueni ri2^* 


84 


d(Cis)/<lt> Limit 


Rate of change in the correction to effluent CO2. 


101 


d(Z pue,)/at > Lmut 


Rate of change in the computed L-Factor. 


103 


d(HHVAF)/dt> Limit 


Rate of change in the computed heating value. 


109 


d(my^)/dt > Limit 


Rate of change in the computed fuel flow. 


111 


d(^H2oydt> Limit 


Rate of change in the computed fuel water fraction. 


121 


d[C2sJAc/(C7sGAct)]/dt 
> Limit 


Rate of change in the parameter of corrected effluent 
H2O divided by the corrected effluent O2. 


131 


d[mAF/HHVAF]/dt > Limit 


Rate of change in the parameter of computed fuel flow 
divided by computed fuel heating value. 



[0080] It is an important aspect of the present invention that it may be integrally involved 

with any of the Input/Loss methods which compute fiiel chemistry. As seen with the used of 
TABLE lA, without a determination of fuel chemistry, i.e., computing quantities based on 
consistent stoichiometrics (with or without tube leakage), then the use of Eq.(200) as taught herein 
to detect tube leaks would become limited. More specifically, this invention is integrally involved 
with The Input/Loss Method of '994 and '879 as tfie determination of fuel chemistry is then based 
on a selection of Choice Operating Parameters which might well require corrections for 
stoichiometric consistencies. For example, if effluent water was being measured, but whose signal 
was not corrected for stoichiometric consistency as taught in '879, resolution of tube leaks (even 
with computed ttj quantities) would be hampered; especially so if fuel chemistry was assumed 
constant. 
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Minimization Techniques. Background 

[0081] The following four paragraphs teach the use of multidimensional minimization 

techniques, including the use of System Effect Parameters and Choice Operating Parameters. 
Taught is how their use allows resolution of flow rates associated with tube failures. Although 
only the Simulated Annealing technique is the Preferred Embodiment for determining tube failure 
flow rate, discussion of all techniques, including neural networks, is made to present the fiill scope 
and spirit of this invention. 

[0082] Any number of multidimensional niinimization techniques may be used by this 
invention; to provide foil teaching four are discussed in detail. All techniques seek to minimize the 
numerical value of an objective fonction. The techniques discussed include: Broyden-Fletcher- 
Goldfarb-Shanno (BFGS), generic Conjugate Gradient, Newton-Raphson and Simulated Annealing 
techniques; references cited below. These techniques, and, notably, their combinations, are 
designed to address all situations of bias in Choice Operating Parameters. All of these techniques, 
except Simulated Annealing, employ derivatives of the objective fonction with respect to the 
independent variable. These techniques all require input of initial values of Choice Operating 
Parameters (AQ.i). The BFGS, generic Conjugate Gradient and Newton-Raphson techniques 
employ unconstrained searches towards optima. Simulated Annealing employs a random but 
constrained search through which tiie Choice Operating Parameters are numerically bounded by 
lower and upper limits. From research and study conducted to develop this invention, the objective 
fonctions described below have proven to be superior for a wide variety of thermal systems burning 
any fossil foels. 

[0083] A common problem facing minimization techniques is the so-called shallow valley 

problem in which an appreciable change in an independent variable has a small effect on the 
objective fonction, even through that change is both real and ^propriate to the physical system. 
This is especially true when applied to the deteraiination of tube leakage flow rate in which a single 
and small tube failure of, say, 5000 Ibm/hr in a large steam generator might represent «0.1% or 
less, of its feedwater flow. Study conducted for the development of this invention, and considered 
unique to it, has found that the Bessel fonction of the first kind is ideally suited to diminish the 
impact of the shallow valley problem. The Bessel fonction emulates the sensitivity that important 
Choice Operating Parameters have on both System Effect Parameters and on the descriptive 
thermal system in general. The Bessel fonction of the first kind of order zero (Jq) has a relatively 
flat (shallow) fimctionality as its argument approaches zero. Apart fix)m this situation, the fonction 
offers non-linearity which is advantageous in converging out-lying arguments. Of great 
importance to those techniques employing derivatives, is that the derivative of Jq is a Bessel 
fimction of the first kind of order one (/j), having a high degree of sensitivity as its argument 
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approaches zero. This derivative relationship addresses a significant number of shallow^ valley 
problems presented by Choice Operating Parameters associated with fossil-fired systems. Another 
technique addressing the shallow valley problem and involving use of the Bessel function is the 
formulation of its argument, termed either A^, X^, X^ or kj [i.e., Joi^i), ^o(^w)» •^(^H) •^(^T)]5 
these arguments are fiiUy discussed below, being defined by Eqs.(202A), (202B), (202C) and 
(202D). 

[0084] The objective ftinction, F, is a fimction of independent variables 3c, i.e., F{x). Of 

uniqueness to this invention, to address inter-dependencies of the Choice Operating Parameters, Xj 
is defined as a scaled Choice Operating Parameter (A^) using the scaling factor Sj; where, initially: 
Xq.! = S'iAq,i ; further discussed in Paragraph 0102. By design, Choice Operating Parameters which 
directly effect system stoichiometrics are used by The Input/Loss Method to compute certain 
parameters which impact the determination of system efficiency. These certain parameters are 
termed System Effect Parameters and, for the Preferred Embodiment, include four general types 
and their associated reference values: the L Factor concept (L^i); the As-Fired fuel flow (m^); the 
higher heating value (HHVj^); and the As-Fired fiiel water fi:action (J^h2o)- Two L Factors are 
presented, L^^^i and Lqq2 as defined by Eqs.(272) & (273). The higher heating value is chosen as 
either: an uncorrected As-Fired value, HHV^; a Dry value, HHYdryJ and/or a MAF value, 
HHVj^. As the Preferred Embodiment, the As-Fired fuel water fiaction is selected only ^en 
determining tube leakage flow rates. However, the system operator may select fi'om any one or 
more or all of these System Effect Parameters (including any one or more or all of the heating 
values), whose differences with respect to reference values are minimized by altering the selected 
Choice Operating Parameters through minimization techniques. The minimization techniques are 
structured to minimize differences between a System Effect Parameter and its corresponding 
"Reference System Effect Parameter^' (teraied: I Vuei-Ref^ ^C02-Ref» "^af-plt^ ^^^1^.^^!^ and 
>FFji2o.Ref)- System Effect Parameters are chosen such that they reflect influences on system 
efficiency tiirough Choice Operating Parameters, and, at the same time, reflect inter-dependencies 
of the Choice Operating Parameters. For example: changes in the concentration of effluent CO2 
(defined as A|s or Aj^) if caused by changes in the fiiel will effect computed fiiel chemistiy, 
computed heating value and computed boiler efficiency, all of which impact system eflSciency. 
However a change in CO2 (measured as a wet concentration) may be caused by a change in the 
concentration of effluent H2O (defined as A2S or A2b)» or a change in combustion air effecting 
effluent O2 (defined as Aj^ or Aj^, or a change in the fiiel's Sodium to Carbon Ratio (defined as 
A9), where any of these changes themselves may directly effect computed fiiel flow and fiiel 
chemistry; thus the importance of inter-dependency of Choice Operating Parameters. 



-47- 

[0085] The following summarizes the objective functionalities for the Preferred 

Embodiment, demonstrating the aforementioned principles: 

Fix) = Eiei/[Si. MKX M^u% ^W] 

Ah = /[HHVk3, HHVk3.R,^ Mh] 

The symbol Y^iei is defined following Eq.(203). Note that as F (jc) is minimized the quantities A 
are updated in tum (Aj = x/Sj), thus allowing System Effect Parameters to be computed leading 
directly to the computation of A^, A^, A^ and Aj. System Effect Parameters have general 
dependency on Reference Fuel Characteristics, including the following important inter- 
relationships: computed fiiel chemistry is dependent on several or all Choice Operating Parameters, 
A; computed heating values (HHV,^, LHVP and HHVP) are dependent on fiiel chemistry, thus A; 
and boiler eflBciency (rig) determined using the methods herein and those of '429, is dependent 
directly on effluents CO2 and O2, is also dependent on fuel chemistry, and is also dependent on 
heating value, thus K. All of these quantities (fiiel chemistry, heating values and boiler efficiency) 
are also dependent on in-leakage of working fluid into the combustion path, the terms b^ and b^. 
Working fluid energy flow (BBTC), the Firing Correction HBC and the Aenthalpy term associated 
with tube failure enei^ flow (m^Ah), are all dependent on Operating Parameters. 

= SjAi 
Lj^i = /[fuel chemistry(A)] 

niAF = /[(BBTC-mjAh), rig (A), HHVP(A), HBC] 
HHVj^ = /[fuel chemistry(A)] 

WFy^Q = /[tube leakage flow rate (Ag), fiiel chemistty(A)]. 
System Effect Parameters 

[0086] As discussed. System Eflfect Parameters include four general types of parameters 

and theh- associated reference values: the L Factor (L^i); the As-Fired fiiel flow (mj^); the higher 
heating value (HHV,^); and the As-Fired fiiel water fraction (WF^q). The most important of these 
is the L Factor, used routinely for most situations. The higher heating value may be employed, for 
example, when the thermal system is operating under controlled conditions (e.g., under a testing 
program), in which its fiiel is well characterized. Also, during initial installation of a Calculational 
Engine, heating value may be used for scoping the range of reasonable correction factors. Fuel 
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flow is discussed below. The As-Fired fiiel water fraction is typically used for the special case of 
computing the tube leakage mass flow rate, m-p. 

[0087] The L Factor is important in reducing the impact of the shallow valley problem 

found with fossil-fired systems. An important reason for this is that both Z pugj and Lqq2 have been 
demonstrated to have remarkably small standard deviations for fossil fuels, including black liquor 
fiiels. For example, some Ranks of coal have L^uqi standard deviations as low as ±0.05%. 
Although black liquor dry chemistries may vary widely depending on the particular process, a 
given process will typically produce a highly consistent L Factor (either I p^^j or Icoi)- Use of 
^ Fuel ^ ^ System Effect Parameter is the Preferred Embodiment, that is when V^^i is computed 
using Eq.(272) or its equivalence. For black liquor fiiels the methods taught in '563 and '879 
require modification to address the fiiel's inorganic compounds, although the principles taught in 
'563 and '879 remain and are applicable. Specifically the L Factor is computed based on 
theoretical combustion of dry fiiel, thus: (p^ = a^^_2 ^ ^dry-2 "^^a^^z^^^ ^Aa ^ ^Act ^ "^Act 
= kBL = v = 0.0;andaN = /?Act=l-0; and (j)^^ = 4>Ref • 

^Fuel - [^DRY-theorNoRY-Fuel ^DRY-theor 4>Ref ) ^Air 

"(JtheorNH20 ^DRY-theor ^MAF-IO ^Ash 

%T|yfcN|sjaCl %(TNO + kBL)NNa2C03 

^NTNANNa2S04 ^T^oNj^coa)] 
/(XtheorNpuelUHVAF) (272) 

In addition to J^Vuei' ^coi factor has shown promise when used as a System Effect Parameter as 
based exclusively on CO2, given in common units of measure as (Ibm diy CO2 efiEluent)/niillion- 
BtUpugi. Note that is the effluent COj based on theoretical combustion of dried fiiel 
(computed with the same philosophy as used to determine Jj^eor)- 

Lqq2 =10^ PtheorNc02 ] / (XDRY-theorNDRY-FuelHHVDRY) (273) 

The following identities have been found usefiil in determining the L Factor, and for reducing the 
solution problem associated with Eq.(19BL): 



^MAF-Aear NmaF-FucI HH^MAF 



" XtheorNpuciKHVAF 

^ ^DRY-theor^DRY-FuelHHVDRY 



(274) 
(275) 
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When black liquor is burned in a bomb calorimeter for determination of its heating value, ideal 
solid products - as assumed when using Eq.(272) or (273) based on Eq.(19BL), consist of the 
following inorganics: Na2C03, Na2S04, NaCl and K2CO3. However, other ideal products could be 
produced depending on the nature of the fuel, its mix of components, combustion conditions and 
the nature of its chemical equilibrium. Whether potassium forms K2CO3, K2SO4 and/or K2S, or 
forms the same compounds as sodium in the same proportions, or whether fiiel potassium is simply 
combined with fiiel sodium as is a common chemistry assumption, does not affect the general 
scope and spirit of the present invention. Three tools may be used to determine such ideal 
compounds as would then be used to modify Eq.(272), die boiler efficiency computations, 
Eq.(19BL), and throughout this disclosure: 1) direct laboratory analysis of the bomb calorimeter's 
residue; 2) operational e?q)erience; and/or 3) modeling ideal combustion of black liquor using the 
software "HSC Chemistry 5" commercially available from Outokumpu Research Oy, Fori, Finland 
(www.HSC@Outokumpu.com). 

[0088] Along with the L Factor, the system operator may also choose, in any combination, 
the plant's indicated fuel flow, the As-Fired heating value, the diy heating value, the MAF heating 
value, and/or the fuel's average water fraction as System EflFect Parameters. Although the system 
operator has complete flexibility, with this flexibility must apply common engineering judgement. 
For example, optimizing effluent water against HHVj^^^ or HHVj^j^y (heating values without 
water) would make little sense given the lack of connectivity. 

[0089] When selecting the system's indicated fuel flow (m^.p^x) as a Reference System 

Effect Parameter, the recovery boiler's operator should proceed with caution. The flow 
measurement of black liquor feeding a recovery boiler requires an assessment as to its consistency 
and whether the flow measurement encompasses the total feed of fiiel. This mvention teaches that 
the minimization techniques may be used to minimize the difference between a computed fuel flow 
(m^) and the system's indicated fuel flow (m;^.pLx), with an off-set Am^^, dirough converged 
Choice Operating Parameters. Thus, the method of this invention allows use of tiie system's 
indicated fiiel flow to aid in the determination of computed fiiel chemistry, fuel heating value and 
system efficiency. Although not required, for many situations it is the Preferred Embodiment that 
use of the system's fiiel flow be accompanied with either L'^^i or Lqq2 of Eqs.(272) or (273), to 
assist with stability and reasonableness of solution. To fijrther enhance stability and reasonableness 
of solution the recovery boiler's operator may option to limit the range of fiiel concentrations 
determined by the methods of this inventioiL The operator may also limit the numerical range of 
the selected Choice Operating Parameters when using the Simulated Annealing technique. Further, 
to aid in addressing a possible error in m^.pLx, both an off-set (Am^^) in this parameter and a 
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Dilution Factor (M^) may be applied to the relationship between and myy^.pL^; see Eq.(202B) 
below. 

[0090] In summary, the process involving the minimization of differences in System Effect 

Parameters, by optimizing Choice Operating Parameters, results in converged Choice Operating 
Parameters Ap.;; correction factors, Cj, are then determined as based on the ratio of the converged 
Choice Operating Parameter to their initial values iA^{). Aq.^ are based on the system's raw 
instrumentation signal, a previous converged solution, an estimate, a guess or as otherwise 
determined. 

Ci - Ap-i/Ao-i (201) 
Mmimization Techniques. Formulations 

[0091] This sub-section presents general discussions of the multidimensional niinimization 

techniques and detailed formulations useful to the plant operator in niinimizing errors in System 
Effect Parameters. 

[0092] The BFGS technique represents a second generation of multidimensional 

minimization techniques. As such, it is considered one of the most robust of techniques for a well 
conditioned problem. The particular BFGS technique employed by The Input^oss Method has a 
superior reputation for convergence. The only input parameters the user need be concerned with 
are the initial relative step-length and the change in the relative step-length. A well-chosen initial 
relative step-length will prevent long iterations (a value of 0.100 to 0,200 is recommended). The 
change in the relative step-length impacts resolution of the shallow valley problem, and may be 
varied until proper convergence patterns are established. A value between 0.010 to 0.040 for the 
change in the relative step-length has been found to be satisfactory when used in conjunction with 
the scaling techniques taught herein. The BFGS technique is the preferred method for use on a 
continuous bases after the problem has been properly conditioned with scaling factors, and 
selections of Choice Operating and System Effect Parameters have been established. These input 
parameters are also applicable to the generic Conjugate Gradient technique. 
[0093] The generic Conjugate Gradient technique represents a first generation of 

multidimensional minimization techniques. For numerical processing reasons the BFGS technique 
has been demonstrated to be superior in to the generic Conjugate Gradient in convergence 
techniques and accuracy. However, there may be situations in which a generic Conjugate Gradient 
may be useful as an alternative once the problem has been conditioned. 

[0094] The Newton-Raphson method is one of the oldest and simplest multidimensional 

minimization techniques. This method requires the objective function's compounded vector 
gradient, resulting in a Jacobian determinant. Generally it will yield an efficient means of 
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convergence but requires reasonable initial Choice Operating Parameters (A^i); however, without 
such reasonableness it may fail wildly. Newton-Raphson is recommended for possible use only 
after the BFGS technique has failed to meet its convergence criteria. It has applicability given its 
use of the Jacobian determinant, through which forming explicit inter-dependencies between 
System Effect Parameters and all Choice Operating Parameters are employed. This assures 
computed dependencies, if such dependencies exist. This intrinsic feature has been found to be of 
importance when resolving certain power plant problems. The Preferred Embodiment is to 
automatically default fix>m BFGS, given failure to meet its convergence (typically due to a lack of 
established inter-dependencies of Choice Operating Parameters) to, first, the Newton-Raphson, and 
then in-tum, given its failure, to Simulated Aimealing. Newton-Raphson may also be used for 
scopmg initial installations of The Input/Loss Method given diflBcult combinations of System 
Effect and Choice Operating Parameters. 

[0095] The Simulated Annealing technique, because it employs a global, constrained search 

methodology, is the Preferred Embodiment for initial study of a new Input/Loss mstallation and for 
determining tube leakage flow rate. It may also be used to assist in the selection of which Choice 
Operating Parameters are best for a particular thermal system. This procedure simulates the 
annealing process of metal, requiring the controlled reduction of a pseudo-temperature (herein 
termed "pseudo-T") to achieve a desired result (i.e., achieving a minimum potential energy of the 
metal's structure when slowly cooled, thus the minimizing of an objective ftmction). This is a 
bmte force approach involving random search; gradients are not used. As a global optimization 
procedure it may move both downhill and upMl (that is, it may move both towards and away from 
local optima), resulting in distinction between different local optima. Conventional optimization 
techniques (BFGS, generic Conjugate Gradient and Newton-Raphson) only move downhill when 
minimizing an objective function. Conventional techniques are blind to a global solution in the 
sense they immediately choose the downhill direction. When addressing fossil-fired combustion 
problems this may lead to optimizing on tiie most sensitive of a given selection of Choice 
Operating Parameters (most likely COj, thus Ajg or Ajg). Distinction between different local 
optima is accomplished by first starting with initial A^,; values, then successively evaluating 
randomly acquired changes, X, but which fall within user-defined step-lengths. Initially this results 
in coarse study of the objective fimction, employing large step-lengths, requiring repeated 
evaluations with seemingly little progress. In the process of choosing X values the algorithm 
generally attempts to move downhill, however it also moves uphill in a probabilistic manner to 
escape local optima. Step-lengths are dynamically chosen such that half of all uphill moves are 
randomly accepted, again helping to ensure that the function escapes local optima. As the 
annealing process proceeds and the algorithm closes on the global optimum, step-lengths decrease 
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as the pseudo-T decreases requiring even more objective function evaluations as the optimum is 
approached. By viewing objective functions in general terais and witfi its ability to move 
probabilistically uphill. Simulated Annealing solves functions that are otherwise difficult to 
resolve, including shallow valley problems associated witii fossil combustion and tube leakages. 
However, with such flexibility comes numerous objective fimction evaluations necessitating long 
computing times. In addition, converged solutions should be re-tested periodically with different 
seeds (i.e., initializations of the random number generator) to assure a global optimum. 
[0096] When applied to fossil-fired combustion, the more sensitive inputs to the Simulated 

Annealing technique include the following: starting point A^.^ values; the number of cycle 
evaluations (5 is recommended); the minimum and maximum values associated with each Aj (i.e., 
defining the region containing the optimum); an initial pseudo-T (0.100 is recommended); and the 
relative change in pseudo-T (i.e., the step-length, 0.010 to 0.020 is recommended). Each of these 
inputs may be established by sensitivity study to assure a robust solution, or as otherwise 
determined. Minimum and maximum Aj values may also be established by review of historical 
system data or through the experience of the system operator. The smaller the range between 
minimum and maximum Aj values, the tighter the search becomes with the final solution becoming 
narrowed. This feature is especially usefiil when As-Fired fiiel flow or As-Fired fuel water fi*action 
are selected as a System Effect Parameters (in combination witii non-unity Dilution Factors, 
and Mj). 

Objective Function and Choice Operating Parameters 

[0097] The following paragraphs present the preferred objective functions and their solution 

methodologies and specify the Choice Operating Parameters which may be employed by 
minimization techniques. As explained, the Bessel fimction is used to define the objective 
function. The Bessel fimction's argument, as taught by this invention, has been chosen to aid in 
addressing the shallow valley problem and in convergence of the minimization techniques. The 
formulations presented produce quantities which may allow numerical inter-dependencies between 
Choice Operating Parameters (A), or not, depending on which minimization technique is chosen. 
This is important for addressing problems in which initial Choice Operating Parameters lie far fi-om 
the optimum. This is also important where more than one System Effect Parameter is chosen 
which may present unique numerical convergence problems. 

[0098] For the BFGS, generic Conjugate Gradient, Newton-Raphson and Simulated 
Aimealing techniques the objective fimction is given by the following. Note that Ml, M^, and 
Mj are real numbers; & Mjj are typically equal to one, while & Mj are typically greater 
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than one. Again, the System Efifect Parameters, L^^, m^y;, HHV^ and WF^^q, are fiinctions of a set 
of Aj. 

\ = [(Acl - ^kl-Ref)/^kl-Ref]''' (202A) 

= [(^AF - mAF-PLT - AmAF)/(mAF.PLT + ^mAF)]^'* (202B) 

= [(HHVk3 - HHVk3.Ref)/HHVk3.Ref]^" (202C) 

^ = [(^^H20 - WFrno-^d/f^^mo-Kcff (202D) 

F(x) = Eiei {Si[1.0-Jo(\)] + Si[1.0-Jo(^w)] 

+ Si [1.0 - Jo(^h)] + Si [1.0 - Jo(At)] } (203) 

In Eq.(202B) Am^p is an ofif-set, or bias, observed in the indicated flow, m^.pLx. In Eq.(203) and 
as used elsewhere, the symbol X i € i indicates a summation on the index i, v^ere i variables are 
contained in the set I defined as the elements of X. For example, assume the user has chosen the 
following: A]^ is to be optimized to mmimize the error in L'^^\ and HHV)^^^, A23 is optimized for 
Ipue, and m^F (Hif/ = 1.40, Am^p = 0.0), A4 is optimized for Z-Fueb ^nd Ayg is optimized for Zpuei- 
Therefore: X = (A,s, Ajs, A4, A-^), I = {Ajs, Ajs, A4, Ayg}, thus x = (x„ xj, X3, x^); Xj = S,A,s; 
X2 = S2A2s; X3 = S3A4; X4 = S4A7b; where Eq.(203) for this example then becomes: 



F(x) = Si{[l.O-JM] + [LO-JM]) 
+ S2{[1.0-/oW] + [l-0-^aw)]} 

+ S3[1.0-Jo(\)] 

+ S4[1.0-Jo(Xl)] 



Derivatives dF/dXi for the BFGS and generic Conjugate Gradient techniques, based on Eq.(203), 
are given by the following: 

8F/axi = aF/(SiaAi) 

+ Ji(Xh) [dXyi /aAj] + 7i(Xr) [dJ^ /dAi] (204) 

where, for example: 

[dX^ /aAj] = Mw [(fflAF - mAF-PLT - AmAF)/(mAF.pLT + ^af)]"^*"^ 
(am^F /aAi) / (mAF-PLT + ^af) ; 
and where Xl, X^, X^^, and ffly^ are taken as average values. Gradients, 6!Fi /obcj , for tiie Newton- 
Raphson method, thus defining the Jacobian determinant, are given by the following: 
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+ Si Ji(Xh) [dx^ /(Sj aAj)] + Si j.iXy) [dXj /(Sj aA^) (205) 

where, for example: 

[dX^ /aAj] = Mw [(ffiAF - m^.pLT) / hiaf-plt l^"^"^ 
(am^F /8Aj) / (niAF-PLT + ^af) • 

[0099] In the Preferred Embodiment, selection of Choice Operating Parameters may be 
made by the system operator from any combination or all of the following: 



= D^^.^; Stack CO2 (with effects from air pre-heater leakage) (2 US) 

AjB = Dp^^ R^^ ; Boiler CO2 (without effects from air pre-heater leakage) (2 UB) 

Stack H2O (with moisture from air pre-heater leakage) (2 1 2S) 
= j R^ct J Boiler H2O (without moisture from air pre-heater leakage) (2 1 2B) 
A3 = AF ; Air/Fuel mass ratio (2 1 3) 

^4 ^ ^Act 5 ^ Pre-Heater Leakage Factor (214) 
^5 ^ ^ Act 5 Concentration of O2 in the boundary air (215) 
^ ^Ls 9 Indicated limestone flow rate (216) 
Ays = GAct = (g-^2oNTNs + aNkBL/2) + ap; 

Stack O2 (with air pre-heater leakage) (2 1 7S) 
Atb = ^Act (g + 2onTns + ONkBL/2); 

Boiler O2 (without air pre-heater leakage) (2 1 7B) 
Ag = m^; Tube leakage flow rate (218) 

A9 = <Pn ; Sodium to Carbon Ratio in the fuel (2 1 9) 
Aio = <Pr; Reduction Efficiency (220) 



The selection of one or more of the Choice Operating Parameters must depend on common 
understanding of recovery boiler stoichiometrics and associated relationships to physical 
equipment Specifically, The Input/Loss Method produces, by employing one or more of the 
minimization techniques (within the ERR-CALC computer program), converged Choice Operating 
Parameters and correction factors C j q)plied to the initial values A^.^. The converged Choice 
Operating Parameters are then used within the Fuel Iterations to produce a computed fiiel 
chemistry, discussed in conjunction with FIG.2 and FIG.3. A monitoring cycle, processing Fuel 
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Iterations, may be scheduled as frequently as desired; each cycle employing correction factors 
produced by ERR-CALC at the same or slower frequency. For example, ERR-CALC could be 
processed (producing updated correction factors) once per day, while Fuel Iterations could be 
processed once every 2 minutes using 15 minute running averages of Operating Parameters. 
[0100] In the above paragraph, the phase "common understanding of recovery boiler 
stoichiometrics and associated relationships to physical equipment" is meant the routine 
knowledge base a system operator should have concerning his/her thermal system. To thoroughly 
teach this invention, examples of such conmion understanding and their associated impacts on this 
invention follow: if limestone (A^) is not used, the system operator would not select limestone flow 
as a Choice Operating Parameter as such a selection would result in an unity correction factor, non- 
convergence, warning messages, and/or a faulted condition produced from ERR-CALC; the 
selection of the Air Pre-Heater Leakage Factor (A4) would not be made if the system uses a tubular 
exchanger which has no air leakage (as designed), and would result in a similar faulted condition; 
the selection of the air/fuel ratio (A3) leading to determination of the fiiel inert fraction, and also 
invoking a constant fuel inert assumption, would not be made as such a selection would result in a 
similar faulted condition; the selection of Boiler CO2 (Ajg), an Air Pre-Heater Leakage Factor 
(A4), and Boiler O2 (Ayg), given that "correcting" the air pre-heater leakage would have no efifect 
on the Boiler-side mix of CO2 and O2, would result in a similar faulted condition; the selection of 
the Sodium to Carbon Ratio in the fiiel (A9) would not be made if periodic fuel analyses indicates a 
constant ratio. 

[0101] The use of the exponents Ml, M^, and in Eqs.(202A), (202B), (202C) and 
(202D), termed Dilution Factors, allows a dilution or dampening of the functionality between 
Reference System EflFect Parameters (L^ij^sf, m^.pLT, HHVj^.Ref and f^H20-Ref) ^nd selected 
Choice Operating Parameters (A). As an important feature of this invention. Dilution Factors allow 
the numerical processes to recognize that Reference System Effect Parameters may themselves 
have bias. Examples of such bias include: Reference Fuel Characteristics having been chosen with 
an out-dated database, biasing the computed reference L Factor; the reference heating value having 
been determined incorrectly, analyzed incorrectly in the laboratory and/or having intrinsic 
uncertainties; and the indicated fiiel flow having serious instrumentation error. However, 
engineering judgement and a valid database may be reasonably anticipated and applied m the cases 
of reference L Factors and reference heating values. Dilution Factors Ml (influencing L^ij^^d 
Mh (influencing HHV|^.Rgf) may be assumed to be imity for most situations as is preferred; or they 
may be based on monitoring experience, sensitivity studies or as otherwise determined. Applying 
engineering judgement and/or a valid database is difficult in the case of the plant's indicated fuel 
flow; which could have bias, thus the use of Am^ in Eq.(202B). Thus (influencing m^F-PLT) 
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should be determined based on results from The Input/Loss Method and the processes of this 
invention, when such results are generically compared to system data. Specifically, M^,^ may be 
adjusted until Input/Loss computed total effluent flow reasonably agrees and/or tracks the 
measured, computed combustion air flow agrees and/or tracks the measured, computed fiiel flow 
agrees and/or tracks the indicated fuel flow, and similar system-wide comparisons. In the context 
of the last sentence, 'tracks" is defined as the computed value trending over time with the 
measured, having a constant off-set. In the case of fiiel water fraction as used directly as a System 
Effect Parameter, WF^q is driven towards a reference value by optimizing tube leakage flow (Ag); 
when so employed. System Effect Parameter WFii20 considered a special case. Application of 
Dilution Factors require that the sense of the bracketed terms of Eqs.(202A), (202B), (202C) and 
(202D) be always positive requiring a reversal of the derivative's sign as appropriate. The Dilution 
Factors Ml, M^, Mjj and Mj are real numbers; M^ & are typically assumed to be unity, while 
M^ & Mj are typically found through sensitivity studies to be non-unity (ranging between 0.90 
and 1.20 for M^, and 1.0 to 1.6 for Mj). 

[0102] In these relationships each Choice Operating Parameter (A|) is scaled with the 
parameter Sj, determined to be suitable for the BFGS, generic Conjugate Gradient and Newton- 
Raphson techniques. Scaling for these methods is important for proper application of this 
invention, as minimization techniques in general are sensitive to variations in the numerical size, 
and imits of measure, of the A^ terms (e.g., for fossil-fired applications, an un-scaled Aj^ may be 
0.14 moles-COj/mole-Dry-Stack-Gas, while an un-scaled A^ may be 22,000 Ibm/hr). It has been 
found that a good initial estimate of Sj may be developed as the inverse of Aj. Further, the 
influence of scaling may be improved by employing a pre-scaling factor, Sj; which may be 
determined as taught in '879, or as otherwise determined by the system operator through sensitivity 
studies. It has been found that Sj for Ag typically of 10,000 works well if Ag units of measure are 
Ibm/hr. However, it has also been found that the Newton-Raphson technique converges quickly 
when optimizing the combination of Ajs or Ajb, and Ag, thus may be used to adjust Sj for Ag until 
appropriate sensitivity is reached between tiie Choice Operating Parameters of effluent water and 
tube failure flow rate (that is when one term does not predominate the other). 



Si ^ Si/Ao.i 
Xi ^ SjAj 



(206) 
(207) 



[0103] When ERR-CALC is executed using either BFGS, generic Conjugate Gradient or 

Newton-Raphson techniques typically 5 to 50 iterations are required for convergence. However, 
when ERR-CALC is executed using Simulated Annealing typically over 1000 iterations are 
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required for convergence. To address the problem of long computing times, associated with any 
niininiization technique, this invention teaches to duplicate within the ERR-CALC program only 
those calculations which effect System Effect Parameters, and to therefore compute System Effect 
Parameters within ERR-CALC (which are then repeated within the Fuel Iterations). This results in 
a considerable reduction in computing time required to evaluate repeated objective fimction 
calculations. Specifically, these duplicated calculations include: principally HEATRATE 
stoichiometrics (which are also used by EX-FOSS); L Factor calculations; heating value 
calculations; and an approximation of the effects changing stoichiometrics and changing heating 
value has on boiler efficiency and thus the effects on computed fiiel flow using Eq.(330A) or 
(330B). In summary, these duplicated calculations determine affects on the System Effect 
Parameters (Lj^i, m^, HHVj^ and WF^q) of a given set of Choice Operating Parameters (A). 
[0104] Applicable references for the preferred minimization techniques include the 
following sources. For the BFGS and the generic Conjugate Gradient techniques the references 
are: D.F. Shanno and K.H. Phua, "Algorithm 500, Minimization of Unconstrained Multivariate 
Functions", ACM Transactions on Mathematical Software, Vol.2, No.l, March 1976, pages 87-94; 
and D.F. Shanno and K.H. Phua, "Remark on Algorithm 500, Minimization of Unconstreuned 
Multivariate Functions", ACM Transactions on Mathematical Software, VoL6, No.2, Decwnber 
1980, pages 618-622. For the Simulating Annealing technique the references are: W.L. Groffe, 
G.D. Ferrier and J. Rogers, "Global Optimization of Statistical Functions with Simulated 
Annealing", Journal of Econometrics, Vol.60, No.1/2, pp.65-100, Jan./Feb. 1994; for its base 
technology see: A. Corana, M. Marchesi, C. Martin and S. Ridella, "Minimizing Multimodal 
Functions of Continuous Variables with the 'Simulated Annealing' Algorithm", ACM 
Transactions on Mathematical Software. Vol.13, No.3, pp.262-280, Sept 1987; for modifications 
to the random number generator RANMAR which is employed by Simulating Annealing see: F. 
James, "A Review of Pseudorandom Number Generators", Computer Phvsics Communications ^ 
V0L6O, pp.329-344, 1990. For the Newton-Raphson technique the reference is: W.H. Press, S.A. 
Teukolsky, WT. Vettering & B.P. Flanneiy, Numerical Recipes in FORTRAN 77> The Art of 
Scientific Computing . Cambridge University Press, Cambridge and New York (1992), Chapter 9.6 
on Newton-Raphson Method for Nonlinear Systems of Equations, and Chapter 9.7 on Globally 
Convergent Methods for Nonlinear Systems of Equations. 

[0105] Additional minimization techniques and teachings of related mathematical 
procedures which may be applied to this invention, are presented in the following: J. Nocedal and 
S.J. Wright, Numerical Optimization^ Springer- Verlag, New York (1999); G.N. Vanderplaats, 
Numerical Optimization Techniques for Engineering Design. McGraw-Hill Book Company, New 
York (1984); and W.H. Press, S.A. Teukolsky, W.T. Vettering & B.P. Flanneiy, Numerical 
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Recipes in FORTRAN 77. The Art of Scientific Computing. Cambridge University Press, 
Cambridge and New York (1992). Other common minimization techniques involving constrained 
or unconstrained searches may also be altematively applied. These include Sequential Linear 
Programming, Direction Set using Powell's method. Simplex method. Downhill Simplex method. 
Simplex method with product form inverse, Quasi-Newton method, and others. Commercial 
products are also available, such as from Lindo Systems, Inc. of Chicago, 111. 
[0106] A further technique applicable to the reduction of errors in Choice Operating 
Parameters lies with use of artificial neural network technology (herem termed ANN). As 
traditionally employed at power plants ANN technology learns through a database how to 
minimize defined parameters: a change in a coal mill's air registration produces an observed result 
in another parameter. Generically, the minimization of defined parameters means to drive such 
parameters in one, but constrained, direction: the lowest combustion air flow (via bias on FD Fans) 
for a given power; the highest boiler eflSciency by minimizing (1 .0 - Tig); etc. However, as used for 
this invention ANN technology is used to correct Choice Operating Parameters such that errors in 
System Effect Parameters are reduced. Explicit determination of fuel chemistry, computed heating 
value, computed boiler eflSciency, consistent mass and eneigy balances, and, in general, explicit 
thermodynamics remain as taught in '994 and '429. ANN may be ^plied to recognize patterns in 
computed System Effect Parameters influenced by causal Choice Operating Parameters. Much like 
the aforementioned (and preferred) techniques, ANN technology may make corrections to initial 
Choice Operating Parameters to achieve a desired result [for example, to minimize the Al, A.^, Ajj 
and/or Xj terms of Eqs.(202A), (202B), (202C) and (202D)]. These corrections, C j, may be based 
on choosing the highest probability a set of Aj will produce the lowest errors in System Effect 
Parameters relative to Reference System Effect Parameters. An advantage to ANN is that such 
corrections are leamed; that is, ANN improves its correlations with an ever increasing database. 
Typically such learning may be done without use of an objective fimction, but not always as in the 
case of object oriented ANN. Specifically, Choice Operating Parameters used to compute fuel 
chemistries and heating values (leading to boiler eflSciency and system eflSciency), may be 
analyzed for their influences on System Effect Parameters, patterns then recognized which would 
lead directly to C j corrections being applied via Eq.(201). Given such corrections. The Input/Loss 
Method would then proceed as described herein, and in '994 and '429 as applicable. 
[0107J Numerous commercial ANN technology software packages are available, for 
example from: NeuralWare of Pittsburg, PA; California Scientific Software of Nevada City, CA; 
The MathWorks, Inc. of Natick, MA; Pegasus Technology of Mentor, OH a subsidiary of KFx, 
Inc; NeuCo, Inc. of Boston, MA; those available from universities; and those to be found on the 
internet. A particularly Eqjplicable ANN technology is available from Computer Associates of 
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Islandia, NY comprising their Neugents technology. In addition, any ANN technology which 
allows for object oriented programming are directly applicable to this invention as such objects, as 
mathematical kernels, may explicitly correct a variety of Choice Operating Parameters employed 
by one of the Input/Loss methods. 

[0108] However, ANN technology is not the Preferred Embodiment given that such 
technology is historically intended for large databases, databases representing processes too 
complex for explicit ttiermodynamics and/or databases those ^^plicable objective functions are 
unknown or otherwise may not be readily discemed. The teachings of the Preferred Embodiment 
of this invention may be applied directly using ANN technologies which have application 
following the general scope and spirit of the present invention. 

Tube Leakage Flow Rate Computations 

[0109] To determine a tube leakage flow rate using the direct application of the 

aforementioned techniques, that is selecting the Choice Operating Parameter for tube leakage flow 
(Ag), in combination with others, as optimized for any or all System Effect Parameters L^i, m^, 
HHV,^, and/or WF^^q has been demonstrated at power plants burning fiiels whose hydrocarbons 
and fuel water are similar to that of black liquor. However the strict application of such techniques 
is not preferred when determining both the flow rate and the location of the leaking tube. The 
preferred technique for determining a tube leakage flow rate is presented in three steps (termed 
"Passes") detailed below. The technique for determining the location of the failed tube within the 
steam generator is detailed by separate sub-section. In summary, for the Preferred Embodiment, 
once a possible Tube Failure Mechanism has been identified (via TABLE lA or TABLE IB): a 
tube leakage flow rate is determined by optimizing Ag, in combination with other Choice 
Operating Parameters except for A|s, A^, A2S and Ajb- Nominal correction factors to effluent 
CO2 and effluent water are obtained fi*om historical evidence. This achieves stoichiometric 
balance, an initial fiiel chemistry and heating value assuming the nominally corrected effluent CO2 
and effluent water are reasonably accurate. Next a final fuel chemistry and heating value are 
determined but this time as influenced by the determined tube leakage flow rate and all routine 
Choice Operating Parameters, except A2S or A2B; that is. Choice Operating Parameters as would be 
routinely selected whose inter-dependencies are now effected by an established tube leakage. As 
defining a "separative analyses process", it was developed to address the situation where effluent 
water, based on either a measurement or an assumption, was being corrected without regard to how 
such a correction might influence other Choice Operating Parameters, especially tube leakage and 
the important effluent CO2. For example, if in correcting a high effluent water signal (whose Aq.2s 
value reflects an actual tube failure) to a lower nominal value, the resultant dry-based effluent CO2 
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may become badly skewed effecting computed heating value. The preferred process first accepts 
the effluent water value using an historically based correction factor, C2s.hist' i ^-, not optimizing on 
Ajs, AjB, or A2B, but optimizing on Ag and all other Choice Operating Parameters. This 
optimization establishes a computed tube leakage flow rate, consistent fiiel chemistry and a heating 
value given a tube leakage. The computed tube leakage could be essentially zero if determined to 
be stoichiometrically consistent. The process then repeats but including A^^ or Ajg and other 
Choice Operating Parameters, again except effluent water, and using the computed tube leakage 
flow rate. This final process then reflects nominal values given the constrained methodology used 
by Simulated Annealing; for example, computed fuel carbon established in-part fi-om Aj^ will lie 
between duAF-A/mn ^MAF-4/max> affluent water is nominally corrected and is consistent with the 
computed tube leakage. 

[0110] The preferred technique, versus a straightforward application of '879, addresses 

several problem areas found during the demonstration of this invention: the marked insensitivity of 
small tube leakages on system stoichiometrics; correction factors being adversely influenced by an 
actual tube leakage, but the resulting effects of converged Choice Operating Parameters on 
stoichiometrics would mask detection of tube leakage; shallow valley problems aggravated by tube 
failures; and marked differences, and possible problems, associated with scaling Choice Operating 
Parameters especially with widely varying tube leakage flow rates (e.g., fi-om 2,000 to 100,000 
Ibm/hr). Further, this separative analyses process lends itself to the determination of the tube leak's 
location within the steam generator. 

[0111] Pass 0 of the Preferred Embodiment represents a typical monitoring cycle using The 

Input/Loss Method, but where the user has optioned for tube failure checking at each execution of 
the ERR-CALC program. A typical monitoring cycle involves minimizing errors in System Effect 
Parameters (Lj^j, m^^, HHVj^ and/or WFi^q) by optimizing a selection of routine Choice 
Operating Parameters s^plicable to the thermal system and its fiiel, but not use of Ag. For such 
routine monitoring m^ = bp = 0.0 is assumed for the stoichiometrics of Eqs.(19BL) and (200). The 
optimization employs any multidimensional minimization technique, and/or ANN technology, 
appropriate to the thermal system and its fuel. This type of monitoring, using routine Choice 
Operating Parameters, is taught in '879 and herein, and typifies normal use of ERR-CALC. This 
process defines Pass 0 logic described in FIG.3. 

[0112] If tube failure checking (the option having been activated in Pass 0) has detected a 

Tube Failure Mechanism (see TABLE lA and TABLE IB), a Pass 1 process is tiien begun to 
determine a computed tube leakage flow rate. Within this Pass 1 routine Choice Operating 
Parameters used in Pass 0 are selected with Ag, but excluding effluent CO2 and effluent water (A13, 
^iB' ^28 ^2b)- Choice Operating Parameter Ag is optimized to drive the As-Fired fuel water 
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fraction, WFyi2o, ^® computed fuel flow, m^, to their respective reference values: ^^H20-Ref 
(^AF-PLT ^af)- Corrections to eflQuent CO2 and effluent water are set to historical values 
associated with established instrumentation experience (or an assumption) not reflective of tube 
failure. For example, effluent Stack water may be equated to C2s-hisA)-2S where C2s-hist ^ 
nominal correction factor and where A^.2s is based on plant data. Corrections to effluent CO2 and 
effluent water are thus held constant for Pass 1. Other correction factors, Cj, established in Pass 0 
are used with initial Aq^i estimates. This process results in a computed tube leakage flow rate 
which satisfies stoichiometric balances of Eqs.(19BL) and (200), providmg a solution in which fiiel 
chemistry and heating value fall within numerical constraints (min/max limits). Note that other 
System Effect Parameters may be used, optimizing Ag, such as HHY^^, however either WF^q or 
m^ have proven to be appropriately sensitive to the stoichiometrics and represent the Preferred 
Embodiment. This process defines Pass 1 logic described in FIG.3. 

[01 13 J After convergence of Pass 1, Pass 2 then re-establishes System Effect Parameters 
with the previous selection of routine Choice Operating Parameters, but excluding effluent water 
(those effects are now replaced by the computed tube leakage flow rate). The Simulating 
Annealing algorithm is preferred; however for fliis Pass 2 BFGS has been observed to be adequate 
if properly scaled. Convergence results in converged Choice Operating Parameters, noting that all 
Choice Operating Parameters which were being used in Pass 0 are now influenced by a tube 
leakage flow rate (m^), including a nominally corrected effluent water. Pass 1 and Pass 2 employ 
Simulated Annealing technique as the Preferred Embodiment given its ability to address: shallow 
valley problems; the possibility of considerably different scaling associated with small to large tube 
failure flows; its constrained search methodology in which the Choice Operating Parameters are 
numerically bound; and that a variety of Choice Operating Parameters which may have wide 
ranging numerical values. This process defines Pass 2 logic described in FIG.3. Note that 
whenever m^ is determined greater than zero, procedures as taught by this invention then 
incorporate the influence of such leakage into combustion stoichiometrics through the bj molar 
quantity determined as follows: 

br = mj (x Naf) / (^mo ^af) (321 A) 

= W (Nh2o mAp) / (x N^p) (321B) 

The bj quantity, through use in Eq.(19BL), then effects boiler eflBciency, computed fiiel flow and 
efficiency computations in the same manner as the b^ quantity as taught in '994 and '429; b^ in 
'994 and '429 being replaced by the quantity (b^ + bj). Thus if the quantities: N^p; the molecular 
weight of the working fluid (Njjjq, water is assumed); m^ (or rtij^^pu; if the computed is not 
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available); the tube leakage moles (by); and resolution of the stoichiometric model of the 
combustion process, then Eq.(321B) may be solved for m^ without use of a minimization or ANN 
technique. 

[0114] To further assist in teaching this invention, TABLE 2 presents a typical scenario of 

routine monitoring, the identification of a possible tube leak, and then the resolution of the tube 
leakage flow rate. In TABLE 2, the second column denotes the selection of Choice Operating and 
System Effect Parameters; for example, "A^g min Vuer ^^^^ that Choice Operating Parameter 
Ajs, see Eq.(21 IS), is selected to minimize the error in System Effect Parameter Vuci- ^® 
"Ais min \ "A2S min I. Vuei" ^nd "A4 min i^Vuei" ^ 0 is typical for the assumed 
thermal system. If measuring Stack O2 instead of Boiler O2 as used in TABLE 2, the selection 
would typically consist of only "A^g min Z^Vuer ''^2S ^Vuei"- 



TABLE 2: Example of Determining Tube Failure Flow Rate using the Preferred Embodiment 



The Thermal System and 
Computational Sequence 


Optimizations 


Comments 


Routine monitoring of a black liquor- 
fired system having high water, low 
& constant inerts; multiple O2 
instruments are used at the Boiler 
having high accuracy; CO2 & H2O 
Stack instruments; uncertain air 
leakage. Defines Pass 0. 


AisminLVuei 
A2sminrFuei 
A4 minLVuei 
A9 minrpue, 


Compute As-Fired fiiel chemistry with 
constant fuel inerts every 3 minutes; 
optimization of Choice Operating 
Parameters every 30 minutes using 
BFGS technique with Tube Failure 
Model option invoked. Historically: 
Cis = 0.96, C2S = 1.02 and C4 = 1.10. 


A possible tube failure has been 
detected given Tube Failure 
Mechanism 51 (see TABLE lA). Use 
historical values for Cj^ & C2s\ 
Aq^ = 500 Ib/hr. Defines Pass 1 . 


A4minrp^^g, 


Compute As-Fired fuel chemistry with 
constant fuel inerts; optimization using 
Simulated Annealing. Results in the 
computed tube leakage flow rate, mj, 
which satisfies stoichiometrics. 


Return to routine monitoring but 
including the computed tube leakage 
flow rate, mj, but excluding effluent 
water (C2s= 1.02). Defines Pass 2. 


AisminlVuei 
A4 minrpue, 


Compute As-Fired fiiel chemistry as in 
Pass 0, resulting in convei^ged Choice 
Operating Parameters associated with 
the computed tube leakage flow rate. 



/// 
/// 
/// 



Tube Leak Location 

[0115] A important feature of this invention is its ability to assess the impact of a tube leak 
on the thermal performance of the system, and where within the steam generator the leak occurs. 
Once a tube leakage flow rate has been determined, its impact on the total eneigy flow to the 
working fluid and on boiler eflSciency may be determined; thus its effects on fiiel flow and system 
efiBcienq^ may be understood If a thermal system's feedwater flow is held essentially constant, 
then a developing tube leak will result in less total energy flow required from the combustion 
gases; i.e., a reduction in the working fluid's energy flow required to meet the same working fluid 
boundary pressures and temperatures. If the working fluid energy flow without tube leakage is 
teraied BBTC, then the actual energy flow, assuming a tube leak, is given by: (BBTC - m^Ah), 
where is the enthalpy difference between the outlet of the last heat exchanger effected by the 
leakage, hLas^. and the first exchanger so affected, h^^^^j (i.e., the heat exchanger in which the leak 
occurs); mjAh is the energy flow lost from the working fluid due to tube leakage. The enthalpy of 
the leaking fluid as it enters the combustion gas path, hg^ggn^j, is assumed, by choice, to be the same 
as the heat exchanger's inlet enthalpy (any location may be applied). When applied to The 
Input/Loss Method of computing boiler eflBciency, the enthalpy of the leaking fluid entering the 
combustion gas path must be properly referenced as taught by Eq.(137) using (h^^^^ " hf-cai)» 
as taught by Eqs.(131) & (132) using (h^^^ - hsteamr)- Quantitative effects on boiler efficiency and 
system heat rate are not obvious and may not be off-setting, they are preferably described through 
use of The Input/Loss Method of '994 and '429; they may also be described through use of other 
Input/Loss methods, or may be described through use of any method which determines boiler 
eflBciency and system heat rate (and therefore system eflBciency) given knowledge of the fuel's 
heating value and other dependent quantities. Computed As-Fired fuel flow and system heat rate 
are then determined by the following, assuming a tube leak (versus no leak associated with 
Eq.(128)): 





= (BBTC - niTAh) / [ t1b.hhv (HHVP + HBC)] 


(330A) 




= (BBTC - mjAh) / [ t1b.lhv (LHVP + HBC)] 


(330B) 


^HHV 


= niAF (HHVP + HBC) /Wo^ 


(331 A) 




= (BBTC - mT-Ah) / (t,b-hhv W^^^p J 


(33 IB) 




= mAF (LHVP + HBC) /W„^ 


(332A) 




= (BBTC - niTAh) / (t,b.lhv 


(332B) 


ISYS-HHV 


= 3412.1416 /iXRHHV 


(333A) 


^SYS-LHV 


= 3412.1416 /^^LHV 


(333B) 
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[0116] It is important to this invention to recognize that the location of the tube failure 

affects the working fluid's energy flow. The typical recovery boiler used in the paper processing 
industry routes the working fluid first through an economizer heat exchanger, then through a series 
of water wall heat exchangers and the boiler drum (if a sub-critical unit). It could happen tiiat 
additional heat exchangers are also employed, routing the working fluid through a primary super- 
heater, through a final super-heater, and, lastly, through a reheat exchanger. Although the 
preceding is typical, numerous configurations of different heat exchanger names and types may be 
found, including, for example: lower economizers, upper economizers, primary secondary super- 
heaters, fmal secondary super-heaters, division walls, fi^ont reheaters, rear reheaters, etc. Such heat 
exchangers need to be identified by order. For example, if a tube leak occurs in an economizer, its 
loss is seen throughout the steam generator (having the greatest impact on working fluid energy 
flow, e.g., Ah = (hReheat-outiet " hsteamT)^ where hsteamT Icss than 490 Btu/lbm. If a tube leak 

occurs in the final reheat exchanger, its loss only affects this last exchanger (having the least impact 
on working fluid energy flow and thus on computed fiiel flow. Ah = (hRgheat-outiet " hsteami)* where 
hsteamT ^^ay be « 1320 Btu/lbm with a hReheatK)utiet 1520 Btu/lbm. This suggests that Eq.(330A) 
or (330B) has an unique solution dependent on the assigned (and actual) location of the tube 
leakage. When using one of the appropriate Input/Loss methods, such dependency on a location of 
the tube leak may be intrinsically a function of computed fiiel chemistry and Firing Corrections, 
and thus, also a fijnction of the resultant heating value and boiler eflBciency. 
[0117] Determination of the location of the tube leak is accomplished by recognizing that 

certain System Effect Parameters are a fimction of tiie working fluid energy flow (as affected by 
tube failure flow rate and its location). The System Effect Parameter of computed fuel flow, m^, 
is a function of (BBTC - m^Ah) through Eq.(330A) or (330B). The System Effect Parameter of 
As-Fired fuel water fitiction, WF^209 is a function of the (BBTC - m^Ah) term through affects on 
boiler eflBciency (t1b-hhv)» heating value (HHV^) and Firing Corrections (HBC). Although not 
obvious. The Input/Loss Method, because it determines fiiel chemistry, heating value, boiler 
efficiency and Firing Corrections independent of fuel flow, and with great consistency, must, 
never-the-less, affect computed boiler eflBciency consistently as taught herein and in *429; 
illustrating using higher heating value: 

^B-HHV = (BBTC - m^Ah) / [m^p (HHVP + HBC)] (334) 
= (- HPR;,^.HHv + HRXa^.hhv) / (HHVP + HBC) (335) 

Eq.(335) must reflect a consistently computed boiler eflBciency; just as Eq.(334) as composed of a 
term which directly reflects tube failure location (BBTC - m^Ah), and therefore reflects a 
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consistently computed boiler efficiency. The Enthalpy of Products and the Enthalpy of Reactants 
terms of Eq.(335), HPR^ct and HRX^ct^ computed with terms influenced by both the tube 
leakage flow and its location via the Ah term. HPRy^^j^ includes the enthalpy of ail water exiting the 
system (hs^ack)? relative to the enthalpy at associated entry points into the combustion gas path, 
including that for tube leaks at hsteg^y; thus: (h^^^ - hstgaj^^) as found in Eqs.(131) & (132). 
HRX^^jt of Eqs.(135) & (136) includes the Firing Correction term which encompasses the entering 
enthalpy of all in-leakages of water, including tube leaks at h^tggjjjj, relative to a reference enthalpy 
taken as the saturated liquid enthalpy at the calorimetric temperature; thus: (hg^gaj^T - hf_Q^ as 
found in Eq.(137). 

[0118] Determination of which heat exchanger has a tube leak is accomplished by assigning 

the tube leak to successive heat exchangers, in repetitive computations involving ERR-CALC and 
Fuel Iterations, and then examining "key comparative parameters" produced from these 
computations for deviations from their "reference key comparative parameter". Reference key 
comparative parameters are determined from a Fuel Iteration witfiout tube leakage, or as otherwise 
obtained. If using the System Effect Parameter of computed fuel flow, m^^, such key comparative 
parameters include: the As-Fired fiiel flow, the average fiiel water fraction, heating value and the 
Fuel Consumption Index associated with that heat exchanger. The following weightings of these 
key comparative parameters has been found useful in determining the lowest deviation among tiie j 
heat exchangers when using System Effect Parameter of fuel flow: 



If using the System Effect Parameter of As-Fired fiiel water fraction, WFu2o, such key comparative 
parameters include: the As-Fired fiiel flow, the average fuel water fraction, heating value and the 
Fuel Consumption Index associated with that heat exchanger. The following weightings of these 
key comparative parameters has been found useful in determining the lowest deviation among the j 
heat exchangers when using System Effect Parameter of fuel water: 



Deviatiouj = 0.02 (m^p - hIaf-plt " ^AF)/(mAF.PLT + ^AF)|j 

+ 0.08 I (JFFH20-Ref - WF^O) I ^H20-Ref I j 

+ 0.05|(HHVRef -HHVAF)/HHVR,f|j 




(382) 



Deviationj = 0.08 ((m^p - mAp-PLT " ^AF)/(mAF.PLT + ^AF)|j 

+ 0.02 I (>^H20.Ref - WF^o) I ^H20.Ref I j 

+ 0.05|(HHVR,f -HHV^)/HHVR,f|: 
-h 0.85|(FCW.FCI)/FClRef|j 



(383) 
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Other key comparative parameters have been studied and have been found as useful as those used 
in Eqs.(382) & (383), depending on the thermal system and its fiiel, and the convergence criteria 
employed by the minimization techniques. These other key comparative parameters include those 
in the following list; however, this list is not meant to be exhaustive but representative of the scope 
and spirit of the present invention: 

Boiler efficiency; 

Dry heating value; 

Specific working fluid energy flow [(BBTC - mjAh) / m^]; 
Weight fraction of fixel carbon as computed by one of the Input/Loss methods; 
Weight fraction of fuel hydrogen as computed by one of the Input/Loss methods; 
Numerical convergence associated with Minimization Techniques including ANN; 
Computed overall heat transfer coefficient of the assumed leaking heat exchanger; 
Computed log-mean-temperature difference of the assumed leaking heat exchanger; 
Computed cleanliness factor of the assumed leaking heat exchanger; and/or 
Computed relative irreversible thermodynamic loss of the assumed 
leaking heat exchanger. 

The Preferred Embodiment of this invention is to employ Eqs.(382) and (383). However, the 
weightings found in these equations are presented to demonstrate the general technique of 
determining a niinimum deviation, and/or otherwise to determine flie location of the failed tube. 
Such weightings may easily change depending on the nature of the fiiel being burned, the design of 
the thermal system, the key comparative parameters employed, and the specific computational 
procedures employed from one of the various Input/Loss methods. These weightings and the 
development of a set of key comparative parameters may be determined through simulations of 
failed tubes, computationally locating the failures within different heat exchangers of the thermal 
system. The concept and development of Fuel Consumption Indices is described in the following 
reference: FD Lang, "Fuel Consumption Index for Proper Monitoring of Power Plants - Revisited", 
Am. Societv of Mech. Engrs.> 2002 International Joint Power Generation Conference, Scottsdale, 
AZ, IJPGC2002-26097. 

[0119] The Preferred Embodiment used to resolve the location of the heat exchanger 

containing a failed tube is to first establish a set of "reference key comparative parameters" 
associated with no tube leakage, then to determine the net energy flow to the working fluid 
assuming the tube leakage is in a particular heat exchanger, repeating such determination until all 
heat exchangers have been so analyzed. Finally, deviations are formed, similar to Eqs.(382) or 
(383), from which the lowest deviation will indicate the heat exchanger with the failed tube. Note 
that if using tiie procedures taught in Paragraphs 0109 through 0114 to determine the tube failure 



-67- 

flow rate, then these same computational methods are used at every evaluation of the net energy 
flow to ttie working fluid and associated assumption of where the tube leakage is located. In 
summary, this invention teaches: to identify a set of heat exchangers descriptive of the thermal 
system as employed to transfer net energy flow to the working fluid from the combustion gases 
resulting in a set of identified heat exchangers; to then obtain a set of Operating Parameters 
applicable to the set of identified heat exchangers; to then determine a set of net energy flows to the 
working fluid from the combustion gases based on the set of identified heat exchangers, the set of 
Operating Parameters and the tube leakage flow, each said set of net energy flows descriptive of 
the thermal system and wherein the tube leakage flow is assigned to a different heat exchanger in 
each said set; to then determine a reference key comparative parameter for the thermal system 
resulting in a type of key comparative parameter; to then obtain a set of key comparative 
parameters associated v^th each identified heat exchanger and based on the set of net energy flows 
and the type of key comparative parameter; to then determine a set of deviations between the set of 
key comparative parameters and the reference key comparative parameter; to then determine an 
identification of the leaking heat exchanger based on the set of deviations; and to then report to the 
operator of the thermal system the identification of the leaking heat exchanger such that corrective 
action may take place. 

Conclusion 

[0120] Although the present invention has been described in considerable detail with regard 
to certain Preferred Embodiments thereof, other embodiments within the scope and spirit of the 
present invention are possible without departing from the general industrial ^plicability of the 
invention. For example, the tube leakage flow rate may be computed by applying the techniques 
discussed in Paragraphs 0076 through 0108 (i.e., using a single computational pass). The tube 
leakage flow rate may be determined by use of any System Effect Parameter which indicates an 
appropriate sensitivity to system effects (with, or in addition to, L^i, m^^, HHV|^ and WF^2o)' The 
separative analyses process described above may acquire numerous variations; for example: Pass 1 
correction factors may be set to unity; or Pass 1 may only employ the Choice Operating Parameter 
of tube leakage flow rate. The method described is applicable to either higher heating value or 
lower heating value computations as supported by the teachings in '429. Further, aforementioned 
descriptions of this invention assume that a steam generator's working fluid is water, however the 
general procedures of this invention may be ^plied to any type of working fluid provided that the 
working fluid is definable at the boundary of the system. Examples of other working fluids are: 
mixtures of water and organic fluids, organic fluids, liquid metals and so forth. Further, the L 
Factor concept as used to optimize Choice Operating Parameters may take numerous forms. 
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although two are demonstrated herein, others are discussed in '563, still others may be formed by 
one skilled in the art based on the teachings herein and in '879 and in '563. Further, the Reduction 
EflBciency parameter, the Sulfur to Smelt ratio and the Sodium to Carbon ratio in fuel are defined 
herein as molar ratios whereas they could be re-defined and incorporated within the teachings 
herein as mass ratios, and/or determined fi-om laboratory analyses or fi-om common mdustrial 
experience as mass ratios and then converted to molar ratios. Accordingly, the general theme and 
scope of the appended claims should not be limited to the descriptions of the Preferred 
Embodiment disclosed herein. 

[0121] Although a Preferred Embodiment of the present invention has been illustrated in 
THE DRAWINGS and described in considerable detail the foregoing DESCRIPTION OF THE 
PREFERRED EMBODIMENT, it will be understood that the invention is not limited to the 
embodiments disclosed, but whose methods are cq>able of numerous rearrangements, 
modifications and substitutions without departing fi-om the scope and spirit of the present invention 
as set forth and defined by the claims herein. 

THE DRAWINGS 

[0122] FIGURE 1 is a schematic representation of a themial system, particularly a recovery 
boiler system illustrating use of stoichiometric relationships important in applying this invention. 
It should be studied in conjunction with combustion stoichiometrics of Eq.(19BL). FIG.l depicts a 
recovery boiler system denoted as 20. In this system 20, a black liquor fiiel feed 22 and 
combustion air 24 are all provided to the upstream side region 26 of the heat exchangers & 
combustion region 28. Note ttiat this region 28 does not include the air pre-heater 36. In addition, 
in some types of recovery boilers 20 such as those using fluidized bed combustors, other materials 
may be injected into region 26, such as a flow of limestone 31 to minimize effluent SO2 by 
chemically binding sulfur as CaS04. Other sorbents may be injected to control sulfijr, to control 
other pollutants, and/or to control the combustion process. The black liquor fuel feed 22 contains, 
in general, combustible organic material, sodium-based inorganic material, water and mineral 
matter (commonly called inerts, or ash). Inerts are an unbumable component that passes through 
the system with little physical change, but which are heated and cooled. In the heat exchangers & 
combustion region 28, the black liquor fiiel 22 is burned with the combustion air 24 to form hot 
products of combustion. Heat from the products of combustion is transferred to a working fluid 
that enters 23 heat exchangers 32 that are depicted as integral with the heat exchangers & 
combustion region 28- The heated working fluid 30 is used in a manner ^propriate to a working 
fluid to generate a usefiil output 33 (for example, in a conventional power plant such usefiil output 
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may be supplied to a turbine cycle for the production of electrical power). Heat exchangers 32 may 
consist of a series of heat exchangers as explained in Paragraph 0116. One heat exchanger of the 
collection of heat exchangers 32 may develop a leakage of its working fluid 27, which leakage 
mixes with the products of combustion 28. There may also be working fluid leakage 29 into the 
products of combustion 28 and into region 35, not associated with water in the fuel feed 22, or heat 
exchanger leakage 27, or moisture in the combustion air 24. Working fluid leakage 29 consists of 
known flows, or flows which may be otherwise reasonably assumed or determined; and may result 
from, for example, soot blowing associated with coal-fired systems, or working fluid used to 
atomize the black liquor fuel 22 before combustion, or used in pollutant control processes located 
at 35 or 42. The smelt 37 from the combustion of black liquor fuel is removed from region 28. 
After smelt removal, the remaining products of combustion leave the heat exchangers & 
combustion region 28 on its downstream region 34, the cooler products of combustion then 
commonly flow through ducts, region 35, which may contain fly ash removal equipment, passing 
then to an air pre-heater 36, where a further portion of the combustion gas energy is transferred to 
an incoming air stream 38, which air then becomes the combustion air 24. The total air delivered 
to 20 is the incoming air flow 25. In many cases, an air leakage flow 40 enters the flow of the 
products of combustion as it passes through the air pre-heater 36. The further cooled products of 
combustion leave the air pre-heater 36 and pass to the Stack 42, then being exhausted to the local 
environment. Withm recovery boiler system 20 the combustion gas path is defined as that region 
encompassing the flow of products of combustion, said products also termed combustion gases, 
generally occupying regions 28, 35, the gas side of 36, and 42. 

[0123] FIG.l, given its general system description provided above, is applicable to a wide 

variety of fossil-fired systems, including recovery boilers, traditional power plants, oil-burning 
power plants, gas-fired power plants, biomass combustors, fluidized bed combustors, conventional 
electric power plants, steam generators, package boilers, combustion turbines, and combustion 
turbines with heat recovery boilers. This list is not meant to be exhaustive, however, and is 
presented to illustrate some of the areas of applicability of the present invention which encompass 
any thermal system burning a fossil fiiel and which has at least one heat exchanger whose working 
fluid is being heated by the products of combustion. This invention is applicable to a wide variety 
of Input/Loss methods, especially its ability to identify the location of the failed tube. If a thermal 
system is to be characterized quantitatively using one of the Input/Loss methods, then relationships 
between Choice Operating Parameters to energy flow inputs and outputs as they are affected by 
tube leakage flow rate 27, may be understood with enhanced accuracy using this invention. This 
understanding, in turn, permits the operation of the fossil-fired system to be optimized for 
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efficiency, reduction of polluting effluents, prevention of black liquor explosions and reduced 
down-time given a tube failure and knowledge of its location. 

[0124] Within fossil-fired systems, some quantities are readily measured with adequate 

accuracy, and others may not be measured on-line (in real time) with accuracy sufficient to 
quantify the operation of the system 20 to the required accuracy to optimize efficiency and/or for 
the detection of a tube failure and its location. For example, working fluid flow rates, pressures 
and temperatures may be readily measured with good accuracy by conventional sensors located at 
defmed boundaries such as 23, 30, 25, 33, 42, 22, 29, 31 and 37. Choice Operating Parameters all 
may, under idea conditions, be directly measured with high accuracy either in real-time or 
periodically. However, if they are not measured with high accuracy, the ability of Input/Loss 
methods to quantitatively improve system efficiency may then be compromised. In FIG.l 
quantities which may be (or are) Choice Operating Parameters include: the combustion gas 
concentrations in the regions 35 and 42 (including CO2, H2O, and O2, termed Ajg, A2B, A73 at 
region 35, and A,s, A2S, Ayg at region 42); the indicated combustion air flow 24 (when combined 
with indicated fuel flow then allows the Air/Fuel ratio to be determined, A3, which allows fiiel inert 
fraction to be computed as taught in '994); the ratio of gas concentrations across the air pre-heater, 
regions 35 and 42 (either the Oj or the CO2 ratio across these regions, preferably the CO2 ratio, 
thus allowing the Air Pre-Heater Leakage Factor i?^^ to be determined, A4); the concentration of 
O2 in the combustion air local to the system 25 (termed A^^ or A5, allowing ^^^^ to be 
determined); the indicated limestone flow 31 (A^); the Sodium to Carbon Ratio in the black liquor 
fuel 22 determined in real-time or periodically (A9); and the Reduction Efficiency term determined 
in real-time or periodically through analysis of the smelt 37 and the black liquor fiiel 22 (Aiq). In 
addition to these nine, another Choice Operating Parameter is tube leakage flow rate 27 (Ag), 
which, in the Preferred Embodiment, is optimized using the fuel's average water content in the fiiel 
(FrFn2o-Ref) or using the computed As-Fired fiiel flow (m^); when optimized, the tube leakage 
flow rate becomes defined consistent with stoichiometrics of Eqs.(19BL) and (200), using in-part 
Eq.(321B). Refer to Eqs.(211S) through (220). This invention teaches how to correct such 
measurements or their assumptions if such measurements are not available; or, in the case of a 
recovery boiler system having a tube leakage, a measurement of such tube leakage flow rate being 
impossible to obtain in any reasonable manner without using the teachings of this invention. 
[0125] FIGURE 2 illustrates an important portion of this invention, specifically the general 

calculational sequences associated with optimizing Choice Operating Parameters and subsequent 
Fuel Iterations when monitoring a recovery boiler system on-line, i.e., in essentially real time. Box 
250 represents general data initialization including establishing Reference Fuel Characteristics, 
data collection, data organization and routine set-ups of all programs. Box 255 depicts the use of 



-71- 

the ERR-CALC program, detailed in FIG.3, which produces converged Choice Operating 
Parameters and, given a heat exchanger leak, the tube leakage flow rate. If ERR-CALC is not to be 
exercised at the same frequency as the Fuel Iterations, Box 255 is bypassed; and, if bypassed, its 
previously computed correction factors are applied to Aq.}, then employed within the Fuel 
Iterations. ERR-CALC may be employed, as-needed, for updating correction factors. Box 260 
depicts the FUEL program which reduces fiiel data from identified multiple sources, prepares a 
composite fuel, and then prepares an input file for the system simulator EX-FOSS. Reduction of 
fuel data involves combining the primary (computed) fiiel chemistry from a previous iteration, with 
secondary fiiels which have constant and known chemistries, producing a composite fuel. Box 270 
is system data supplied to the process as on-line (or real time) including at least the following 
Operating Parameters (refer to Paragraph 0031 for details): working fluid pressures, temperatures 
and flows, air psychrometrics, useful system output and other related data. Box 280 depicts the 
system simulator EX-FOSS which, given specification of a composite fiiel from FUEL, inputs 
from Box 270, routine set-up data and converged Choice Operating Parameters from Box 255 
(including possible tube leakage flow rate), produces the following: boiler efficiency using the 
methods of '429 as modified by the present invention, As-Fired fixel flow (m^) using Eq.(330A) 
or (330B), complete effluent concentrations, system heat rate using Eq.(331A), (33 IB), (332A) or 
(332B)5 effluent flow, emission rates of all effluents including the common pollutants, and other 
thermal performance parameters including, for example, energy flow to the working fluid given a 
tube leak (BBTC - m^-Ah), and the Firing Correction (HBC). The determination of many of these 
parameters is taught herein and in '994. Box 285 depicts the HEATRATE program within which, 
given the converged Choice Operating Parameters (including possible tube leakage flow rate), 
produces fiiel chemistry, L Factors and fuel heating value for both the composite fiiel (as either 
higher or lower heating values), and, given the known compositions of secondary fuels, the 
composition of the primary (unknown) fuel. Box 285 also depicts the computations within 
HEATRATE which lead to identification of which heat exchanger has the tube leakage; refer to 
Paragraphs 0115 through 0119. Designation 287 tests for convergence of the process based on 
composite fiiel moles (x), certain effluents such as CO2 and H2O, heating value and computed fiiel 
water fiction; if the convergence criteria is not met the process continues to iterate fi-om Box 260. 
In general, convergences lie within 0.5x10"^ percent of the computed As-Fired fiiel moles. Note 
that the iterations encompassing 260, 270, 280, 285 and 287 define what is meant by the term "Fuel 
Iterations". In summary, fuel Iterations are the iterative calculations between EX-FOSS, as input 
with known fuel chemistry and heating value fi-om a previous iteration but with imknown effluents 
(to be computed by EX-FOSS, except for effluent O2 which is input), and HEATRATE as input 
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with known effluents (i.e., the corrected Choice Operating Parameters) but with unknown fuel 
chemistry and heating value (to be computed by HEATRATE). 

[0126] Designation 290 of FIG.2 illustrates an important teaching of this invention as it 

represents a decision to either: 1) continue towards reporting results and quitting if either the Tube 
Failure Model (i.e., this invention) is not invoked, or is invoked and all computations and processes 
have been completed; or 2) if the Tube Failure Model is invoked and the location of the failure is 
desired, to then continue the computations, starting afresh with Box 255 but with every pass 
through Box 255 to assume the tube leakage occurs in a different heat exchanger. If the tube 
failure flow rate and its location are both desired, then summary analysis of these computations 
(leading to identification of the leakage's location) is performed wdthin the HEATRATE program, 
Box 285, after all heat exchangers have been analyzed. Typically, monitoring cycles are scheduled 
for every 2 to 3 minutes using updated data based on 15 minute running averages. Once conveiged 
and all computations have been completed. Box 294 produces reportable results from the EX- 
FOSS and HEATRATE programs, including system heat rate, tube leakage flow and its location, 
and other thermal performance parameters which are influenced by the failed tube including: 
Second Law analysis of the thermal system (producing Fuel Consumption Indices), fuel flow, total 
effluent flow, emission rates, other output and reports to system operators as to what corrective 
actions may take place; smd reports also being provided to regulatory authorities as requested or 
required. Box 298 of FIG.2 is to quit 

[0127] FIGURE 3 illustrates another important portion of this invention, specifically the 

organization of the ERR-CALC program used to determine correction factors to the initial Choice 
Operating Parameters, including the computation of tube leakage flow rate if warranted. In FIG.3 
Box 310 depicts the start of the program which invokes data collection and routine program set-up 
associated with ERR-CALC. Routine program set-up associated with ERR-CALC includes a user 
option as to whether the Tube Failure Model will be invoked, or not, for a given monitoring cycle 
(also see Box 250 of FIG.2). Box 320 depicts initializations of data including organization of data 
arrays associated with selected Choice Operating and System Effect Parameters, and determination 
of scaling factors, Sj, and pre-scaling factors, Sj. Box 330 depicts the selection, collection and 
processing of general input data associated with the minimization techniques, principally the 
selection of Choice Operating Parameters, System Effects Parameters, Reference Systems Effects 
Parameters and routine inputs, options and convergence criteria to the minimization techniques as 
are known to those skilled in the art using these techniques (such inputs and criteria are presented 
in cited references in Paragraphs 0104, 0105, and 0107). 

[0128] Within FIG.3, if the Tube Failure Model is optioned. Boxes 340, 350 and 360 may 

be executed three times for the Preferred Embodiment; these are termed Passes 0, 1 and 2. The 
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initial execution of Boxes 340, 350 and 360 as a series is termed Pass 0. Box 340 depicts 
application of the minimization techniques as herein discussed, including evaluation of an objective 
function resulting in optimizing the selected Choice Operating Parameters. Box 350 depicts the 
use of a simulation principally of the HEATRATE program within ERR-CALC by which the 
computing time required for the supporting computations required for Box 340 are greatly reduced; 
refer to Paragraph 0103 for details. For Simulated Annealing, ANN and other such exhaustive 
procedures, having importance in implementing the Preferred Embodiment of this invention given 
a tube failure. Box 350 is typically caused to be executed from Box 340 thousands of times for 
each Pass, hiputs to Box 350 are principally Choice Operating Parameters. Output from Box 350 
to Box 340 being principally System Effect Parameters from which the objective function is flien 
evaluated. After convergence of the minimization techniques of Box 340 and thus resolution of 
converged Choice Operating Parameters, Box 360 depicts the calculation of correction factors 
associated with the selected Choice Operating Parameters. Box 360 also includes the production of 
appropriate warning messages associated with the ERR-CALC computations; for example: non- 
convergence, computational failures, the automatic switching to altemative minimization 
techniques, and the like. If the monitoring cycle is processing Pass 0 and the Tube F^lure Model 
is not optioned, the program quits via Box 470. 

[0129] If the Tube Failure Model is optioned, logic flows as discussed above, from Box 340 

through Box 360 resulting in a converged solution (but without an assigned tube leakage). 
However, within Box 360 logic is then invoked which initiates Pass 1 or Pass 2, either proceeding 
to Box 430. For Pass 1, Box 430 tests for a trip mechanism of TABLE 1 A and TABLE IB, and if 
found that a tube failure is possible, then: re-sets A^^ or Ajg, and Ajs or A2B correction factors to 
historically based factors (or to unity); selects the System Effect Parameter for fuel water (^h2o) 
or for fuel flow (m^^) whose difference with reference values is reduced by optimizing Choice 
Operating Parameter for tube leakage flow (Ag); selects other Choice Operating Parameters and 
associated System Effect Parameters, excluding Ajg, A13, A2S and A2B; selects Simulated 
Annealing as the minimization technique; employs the computed tube leakage flow (m-j-); and 
passes control to Box 340. Convergence via Boxes 340 and 350 results in a stoichiometric 
consistency with the computed tube leakage flow rate. After Pass 1, Box 430 logic then initiates 
Pass 2 proceeding to Box 430. For Pass 2, Box 430 then: re-sets all correction factors to unity; 
selects the same System Effect Parameters and Choice Operating Parameters established in Box 
330 (excluding A2S and Ajb), plus System Effect Parameter for fuel water or fuel flow whose 
differences with their reference values is reduced through optimization; selects the Simulated 
Aimealing technique; employs the computed tube leakage flow rate of Pass 1; and passes control to 
Box 340. Convergence via Boxes 340 and 350 results in correction factors for all selected Choice 
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Operating Parameters which, in combination with the determined tube leakage flow rate via 
Eq.(321BX achieves stoichiometric consistency as stated by Eqs.(19BL) and (200). Within Box 
360, logic produces appropriate warning messages that a tube failure has occurred (assuming the 
tube leakage flow rate is greater than zero), its mass flow rate, the associated Tube Failure 
Mechanism, correction factors, and other routine messages. The execution of ERR-CALC is then 
terminated with Box 470. 

[0130] FIGURE 4 illustrates a plot of both actual and computed water flow entering a 
combustion space of a large commercial steam generator under test conditions. The actual flow is 
measured soot blowing flow which was declared an unknown to the Calculational Engine 
operating The Input/Loss Method as disclosed herein (b^ = 0.0). This test emulated a tube failure. 
The observed sensitivity in FIG.4 is less than 0.1% of feedwater flow thus demonstrating the 
practicality of this invention when applied to recovery boilers. FIG.4 is an outstanding example of 
the sensitivity and accuracy of this invention, demonstrating many of this invention's objectives. 
[0131] Symbols within equations may have been italicized pursuant to Patent OflBce 

publication practices. As used in FIG.l, FIG.2 and FIG.3, and throughout the DESCRIPTION OF 
THE PREFERRED EMBODIMENT, mathematical symbols typed in italics and mathematical 
symbols typed in non-italics have the same meaning when taken in context. Understanding context 
is afforded through the use of subscripts and/or through the normal flow of mathematical 
development. Thus, for example, the symbols: ^^ct* ^Act^ Cp, /ij, Jq, Jj, 

•4iech' ^Vuei^ ^C02j ^af^ defined in Paragraphs 0027 through 0030 and 

elsewhere herein, have the same meaning as, respectively: A^^f ^acv f(X Cj, Cp, hj, h, 

^b^Mech' L'pueb ^coi^ ^AF? AF, HR. The use of italic symbols is used for writing 

style. For example, the lower-case / (or f) describes a functional relationship, defined in Paragraph 
0028, and as used in Paragraphs 0033, 0059, 0064, 0069, 0077 and 0085. For example, as defined 
in Paragr25)h 0027, the letter J (or J) with the subscript "Act" describes the "Total eflQuent water at 
the system's boundary 0 + b^^P); moles/base"; as opposed to as the letter J (or J) with the subscript 
"0" or "1", defined in Paragr^h 0028, describing the Bessel fimction of the first kind of order zero 
or one, and used in Paragraphs 0083, 0085 and 0098; as opposed to J^ech which is defined in 
Paragraph 0075 as the mechanical equivalent of heat. For example, the letter "x" (and related 
terms x^^, ^DRY-theor ^MAF-theor) ^o describe the moles of fiiel per base, defined in 

Paragraph 0027; whereas the terms and x are used to describe generic independent variables 
associated with multidimensional mininiization techniques, defined in Paragraph 0028. The moles 
of fiiel is used principally throughout Paragrs^hs 0044 through 0075; while the terms jcj and x are 
exclusively used in Paragraphs 0081 through 0085, 0098 and 0102. For example, C j is used as a 
correction factor to Choice Operating Parameters (e.g., C2s), defined in paragr^h 0028 and used in 
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paragraphs 0079, 0090, 0099, 0106, 0109, 0112 and 0114; whereas Cp is heat capacity defined in 
paragraph 0029 and used in paragraphs 0071 and 0073. Note that correlation constants denoted, 
for example as A3, B15, B3, C5, etc. and used for the mathematical development of Eqs.(42), (43), 
(48), (49) and (50) are discussed in paragraphs 0056, 0057, 0058, 0077 and 0078. For example, the 
term "h" is the moles of effluent nitrogen without air leakage per base as defined in Eq.(19BL) and 
exclusively used in system stoichiometrics of paragraphs 0044 through 0062; whereas h j is a 
specific enthalpy term defined in paragraph 0029 used in energy equations following paragraph 
0063. 



